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1  Introduction 


1.1  Background 

Multi- anchored  or  tieback  wall  systems  are  often  used  for  temporary 
support  of  excavations  that  have  space  restrictions  due  to  adjacent  struc¬ 
tures,  highways,  railroads,  etc.  In  some  cases,  multi-anchored  systems  may 
remain  as  permanent  structures  after  construction.  In  Corps  of  Engineers 
projects,  permanent  tieback  wall  systems  are  used  as  guide  walls  and 
approach  walls  on  navigation  projects,  and  as  retaining  walls  on  highway 
and  railroad  protection  and  relocation  projects. 

The  behavior  of  multi-anchored  systems  may  be  strongly  influenced  by 
factors  such  as  the  sequence  of  excavation  and  installation  of  anchors,  and 
by  fluctuations  in  the  water  table.  Therefore,  to  obtain  accurate  predictions 
of  the  magnitudes  of  stresses  and  deformations  in  the  structure  and  the  sur¬ 
rounding  soil,  it  is  necessary  to  perform  soil- structure  interaction  (SSI) 
analyses  that  model  the  construction  and  operation  stages  of  the  system. 

For  such  analyses,  adequate  models  for  soils  and  soil-to-structure  inter¬ 
faces  are  required. 

A  substantial  amount  of  research  has  been  performed  in  recent  years  on 
lock  walls,  which  are  an  important  type  of  earth-retaining  structure  for  navi¬ 
gation  projects.  These  studies  have  included  SSI  analyses  of  the  Red  River 
Lock  and  Dam  No.  1  (Ebeling  et  al.  1993;  Ebeling  and  Mosher  1996;  and 
Ebeling,  Peters,  and  Mosher  1997),  the  North  Lock  Wall  at  McAlpine 
Locks  (Ebeling  and  Wahl  1997),  and  Locks  27  (Ebeling,  Pace,  and  Morri¬ 
son  1997),  and  are  good  examples  of  available  state-of-the-art  techniques. 
These  studies  showed  that  the  behavior  of  the  soil-structure  interface  has  a 
significant  influence  on  the  magnitudes  of  the  loads  acting  against  lock 
walls.  They  also  illustrated  that  the  preconstruction  and  postconstruction 
stress  paths  followed  by  interface  elements  are  complex,  often  involving 
simultaneous  changes  in  normal  and  shear  stresses,  as  well  as  unloading¬ 
reloading  due  to  postconstruction  fluctuations  of  the  groundwater  level. 

Gomez,  Filz,  and  Ebeling  (2000)  developed  an  extended  hyperbolic 
model  for  interfaces  and  implemented  it  into  the  finite  element  program 
SOILSTRUCT-ALPHA.  The  model  is  based  on  the  Clough  and  Duncan 
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(1971)  hyperbolic  formulation,  which  was  extended  to  model  a  variety  of 
stress  paths.  Gomez,  Filz,  and  Ebeling  (2000)  performed  a  series  of  inter¬ 
face  tests  between  uniform,  fine  sands  and  concrete.  Some  of  these  tests 
followed  complex  stress  paths  that  included  unloading-reloading  and  simul¬ 
taneous  changes  in  normal  and  shear  stresses.  They  also  earned  out  a  pilot- 
scale  lock  wall  simulation  that  modeled  placement  and  compaction  of  the 
backfill,  surcharge  application,  and  changes  in  the  elevation  of  the  water 
table  behind  the  wall.  By  comparing  model  predictions  to  interface  test 
results,  and  results  of  SOILSTRUCT- ALPHA  analyses  to  measurements 
from  the  lock  wall  simulation,  they  concluded  that  the  extended  load / 
unload/reload  hyperbolic  model  may  provide  accurate  estimates  of  the 
response  of  backfill-to-lock  wall  interfaces. 

Important  similarities  exist  between  the  types  of  loading  that  occur  at 
structure-to-soil  interfaces  in  both  multi-anchored  systems  and  lock  walls. 
Therefore,  it  is  possible  that  the  Gomez-Filz-Ebeling  model  for  interfaces 
could  also  be  used  for  SSI  analyses  of  multi-anchored  systems.  However, 
the  model  was  developed  based  on  the  results  of  interface  tests  performed 
using  uniform  fine  sands,  and  the  model  performance  has  not  been  evalu¬ 
ated  using  coarser  soils. 

The  purpose  of  this  investigation  is  to  evaluate  the  accuracy  of  the 
extended  hyperbolic  model  in  predicting  the  response  of  the  interface 
between  concrete  and  coarse  sand.  A  series  of  virgin  shear  tests  were  per¬ 
formed  under  constant  stress  at  the  interface  between  a  coarse  sand  and  con¬ 
crete.  The  results  of  these  tests  were  used  to  determine  the  hyperbolic 
parameter  values  of  the  interface  following  the  recommendations  given  by 
Gomez,  Filz,  and  Ebeling  (2000).  An  interface  test  was  performed  follow¬ 
ing  a  complex  stress  path  that  included  unloading-reloading  as  well  as 
simultaneous  changes  in  shear  and  normal  stresses.  The  interface  response 
measured  during  this  test  was  compared  to  the  response  calculated  using 
the  extended  hyperbolic  model.  It  was  found  that  the  Gomez-Filz-Ebeling 
interface  model  provided  accurate  estimates  of  the  response  of  this  type  of 
interface.  Therefore,  it  can  be  concluded  that  the  extended  hyperbolic 
model  can  be  used  for  prediction  of  the  response  of  interfaces  between  con¬ 
crete  and  a  variety  of  granular  soils.  The  hyperbolic  parameter  values  of 
the  interface  tested  also  add  to  the  database  of  interface  properties  avail¬ 
able  in  the  literature.  The  extended  hyperbolic  model,  together  with  the 
interface  data  that  have  been  generated,  provide  useful  tools  for  analyses  of 
multi-anchored  retaining  systems  and  other  U.S.  Army  Corps  of  Engineers 
structures. 


1.2  Common  Types  of  Multi-Anchored  Systems 

Multi-anchored  systems  can  be  constructed  using  different  materials  and 
configurations.  The  following  are  the  most  common  types  found  in  practice: 
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•  Vertical  sheet  pile  systems  with  wales  and  posttensioned  tieback 
anchors. 

•  Soldier  beam  systems  with  wood  or  reinforced  concrete  lagging  and 
posttensioned  tieback  anchors. 

•  Secant  cylinder  pile  systems  with  posttensioned  tieback  anchors. 

•  Continuous  reinforced  concrete  slurry  wall  systems  with  postten¬ 
sioned  tieback  anchors. 

•  Discrete  concrete  slurry  wall  systems  (soldier  beams  with  concrete 
lagging)  with  posttensioned  tieback  anchors. 

Figure  1-1  illustrates  the  use  of  a  multi-anchored  system  for  a  typical 
navigation  project.  Because  of  the  space  restrictions  imposed  by  an  adja¬ 
cent  railroad,  excavation  for  the  expansion  of  the  waterway  requires  the 
use  of  a  multi-anchored  system.  For  simplicity,  it  is  assumed  that  the 
multi-anchored  system  depicted  in  the  figure  corresponds  to  a  continuous, 
reinforced  concrete  slurry  wall  with  tieback  anchors.  Tiebacks  consist  of 
posttensioned  tendons  with  a  grouted  anchor  region.  A  berm  of  granular 
material  or  riprap  is  placed  at  the  toe  of  the  wall  to  minimize  erosion  and 
improve  stability. 

Figure  1-2  illustrates  the  typical  construction  sequence  of  a  reinforced 
concrete  slurry  wall.  Initially,  a  trench  is  excavated  using  a  clamshell-type 
tool.  The  excavation  is  stabilized  by  the  use  of  mud  slurry.  The  finished 
trench  acts  as  formwork  for  the  reinforced  concrete  panel.  Placement  of  the 
concrete  using  a  tremie  pipe  displaces  the  mud  slurry  and  leaves  a  struc¬ 
tural  concrete  wall  that  can  be  excavated  and  tied  back  in  much  the  same 
manner  as  the  other  tieback  wall  systems.  The  walls  are  reinforced  using 
preassembled  cages,  which  are  dropped  into  the  slurry  trench  just  before 
concrete  placement.  Slurry  wall  systems  are  usually  0.6  to  0.9  m  thick  and 
can  be  placed  to  depths  of  30  m  or  more.  The  construction  process  can  be 
summarized  as  follows: 

a.  Guide  walls  are  constructed  to  facilitate  positioning  and  alignment  of 
the  clamshell  during  the  excavation  process.  To  stabilize  the  excava¬ 
tion,  mud  slurry  is  kept  inside  the  excavation  to  a  level  above  the 
water  table.  As  illustrated  in  Figure  l-2a,  the  excavation  for  each 
panel  follows  a  staggered  sequence.  Two  end  excavations  are  per¬ 
formed  first,  leaving  a  central  core  intact.  After  the  end  excavations 
are  completed,  the  central  core  is  removed. 

b.  A  stop  end  tube  is  placed  at  one  end  of  the  panel  excavation.  This 
tube  is  extracted  after  concrete  placement  leaving  a  semicircular 
indentation.  This  indentation  serves  as  a  guide  for  the  excavation  of 
the  adjacent  panel  and  allows  the  creation  of  a  shear  key  between 
the  panels. 

c.  Once  the  panel  has  been  excavated  to  the  desired  depth  and  the  slurry 
cleaned  of  fine  excavation  material  (desanded),  the  reinforcement 
cage  is  lowered  into  the  excavation. 
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.  Use  of  a  multi-anchored  system  in  a  typical  navigation  project 
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Figure  1-2.  Typical  construction  sequence  of  a  reinforced-concrete  slurry  wall  (Maurseth  and  Sedey  1992) 
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d.  One  or  more  tremie  pipes  are  used  to  place  the  concrete  without  con¬ 
tamination  from  the  slurry. 

e.  Once  the  wall  is  finished  and  the  concrete  reaches  its  desired 
strength,  the  excavation  and  tieback  installation  process  can  begin. 

Construction  of  the  second  navigation  lock  at  Bonneville  Lock  and  Dam 
required  the  use  of  concrete  slurry  walls  to  retain  the  foundation  of  an  adja¬ 
cent  railroad  line.  Detailed  descriptions  of  construction  procedures  for  the 
continuous  reinforced-concrete  slurry  wall  and  for  the  discrete  slurry  wall 
systems  used  at  Bonneville  Lock  are  presented  by  Munger,  Jones,  and 
Johnson  (1990,  1992)  and  Maurseth  and  Sedey  (1992),  respectively. 


1.3  Response  of  Soil-to-Wall  Interfaces  in 
Multi-Anchored  Systems 

A  waterways  expansion  project,  such  as  that  presented  in  Figure  1-1, 
requires  performing  SSI  analyses  to  determine  the  magnitude  of  the  defor¬ 
mations  of  the  soil  above  the  excavation,  and  the  bending  moments  and 
stresses  in  the  retaining  wall.  Such  analyses  require  close  modeling  of  the 
construction  stages  of  the  multi-anchored  system,  as  well  as  adequate  con¬ 
stitutive  models  for  the  soil  and  for  the  interfaces  between  soil  and  struc¬ 
tural  components.  The  finite  element  analyses  performed  by  Mosher  and 
Knowles  (1990)  for  the  tieback  walls  at  Bonneville  Lock  and  Dam  are  a 
good  example  of  the  available  techniques  that  can  be  used  in  SSI  analyses 
of  multi-anchored  systems. 

Figure  1-3  illustrates  some  of  the  construction  and  operation  stages  of 
the  hypothetical  navigation  project  shown  in  Figure  1-1.  For  simplicity,  it 
is  assumed  that  construction  is  performed  in  the  dry.  After  completion  of 
the  continuous  reinforced  concrete  slurry  wall  (Figure  l-3b),  the  soil  in 
front  of  the  wall  is  excavated  to  an  elevation  slightly  below  the  position  of 
the  first  row  of  anchors.  The  anchors  are  then  installed  and  tensioned 
according  to  the  project  specifications  (Figure  l-3c).  Once  these  anchors 
are  tensioned  and  tested,  excavation  continues  until  reaching  the  position 
of  the  second  row  of  anchors.  The  process  is  repeated  until  reaching  the 
bottom  of  the  excavation  (Figures  1  -3d  and  l-3e).  Once  the  excavation  is 
completed,  the  granular  toe  berm  is  placed  against  the  toe  of  the  wall 
(Figure  l-3f).  During  operation  of  the  navigation  facility,  the  water  level 
outside  the  wall  reaches  its  normal  elevation,  which  may  fluctuate  peri¬ 
odically  during  the  life  of  the  structure  (Figure  l-3g). 

Figure  1-4  illustrates  the  type  of  loading  expected  to  occur  on  a  soil-to- 
wall  interface  element  during  construction  of  such  multi-anchored  systems. 
Immediately  after  construction  of  the  slurry  wall,  the  interface  element  is 
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a.  Initial  state 


e.  Final  excavation  and  installation  of 


Figure  1-3.  Stages  of  construction  that  can  be  modeled  in  SOILSTRUCT-ALPHA  analyses 
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subjected  to  a  normal  stress  Gn  (Figure  l^a).1  Because  little  or  no  relative 
movement  has  taken  place  between  the  soil  and  the  wall,  the  shear  stress 
acting  on  the  interface  element  at  this  stage  may  be  assumed  zero. 
Excavation  of  the  soil  in  front  of  the  structure  may  induce  outward  defor¬ 
mations  of  the  wall  and  reduction  in  the  lateral  stresses  within  the  soil  mass 
behind  the  wall.  The  relative  settlement  of  the  soil  behind  the  wall,  which 
may  take  place  as  a  consequence  of  the  reduction  in  lateral  stresses, 
induces  shear  on  the  soil-to-wall  interface  element  (Figure  1  -4b).  Sub¬ 
sequent  installation  and  tensioning  of  a  row  of  anchors  may  increase  the 
normal  stresses  acting  on  the  interface  element.  Tensioning  of  anchors  may 
also  induce  relative  heave  of  the  soil  mass  behind  the  wall  and  a  conse¬ 
quent  reduction  in  the  shear  stress  acting  on  the  interface  (Figure  l-4c). 

Subsequent  stages  of  excavation  and  installation  of  anchors  may  pro¬ 
duce  progressive  shearing  of  the  interface  element  under  varying  normal 
stress  with  intermediate  cycles  of  unloading-reloading.  During  the  life  of 
the  structure,  fluctuations  of  the  water  level  on  both  sides  of  the  wall  may 
induce  further  cycles  of  unloading  and  reloading  of  the  wall-soil  interface. 

The  type  of  loading  imposed  on  the  soil-to-wall  interface  of  a 
multi-anchored  system  may  differ  from  the  simplified  loading  mechanisms 
illustrated  in  Figure  1-4.  Factors  such  as  the  sequence  of  excavation,  the 
distribution  of  anchors,  the  stiffness  of  the  wall,  and  the  response  of  the 
foundation  soil  will  influence  the  behavior  of  the  soil-wall  system. 


1 .4  The  Gomez-Filz-Ebeling  Interface  Model 

A  number  of  interface  constitutive  models  have  been  developed  by  dif¬ 
ferent  authors.  Quasi-linear  models  have  been  used  by  Goodman,  Taylor, 
and  Brekke  (1968);  Desai,  Muqtadir,  and  Scheele  (1986);  Matsui  and  San 
(1989);  and  Wong,  Kulhawy,  and  Ingraffea  (1989).  Nonlinear  models  have 
been  used  by  Clough  and  Duncan  (1971);  Zaman,  Desai,  and  Drumm 
(1984);  and  Desai,  Drumm  and  Zaman  (1985),  among  others. 

Clough  and  Duncan  (1971)  developed  the  hyperbolic  model  for  inter¬ 
faces.  This  model  has  been  used  extensively  in  SSI  analyses  and  design  of 
geotechnical  structures,  including  analyses  of  lock  wall  behavior  (Ebeling 
et  al.  1993;  Ebeling  and  Mosher  1996;  Ebeling,  Peters  and  Mosher  1997; 
Ebeling  and  Wahl  1997;  and  Ebeling,  Pace,  and  Morrison  1997).  The  hyper¬ 
bolic  model  can  provide  an  accurate  approximation  of  the  interface 
response  under  monotonic  loading  at  constant  normal  stress.  A  recent  study 
by  Gomez,  Filz,  and  Ebeling  (2000)  has  shown  that  the  original  Clough  and 
Duncan  (1971)  interface  model  is  not  accurate  for  modeling  the  interface 
response  under  simultaneous  changes  in  shear  and  normal  stresses  or  for 
unloading-reloading . 


1  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  C). 
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Gomez,  Filz,  and  Ebeling  (2000)  introduced  the  extended  hyperbolic 
model  for  interfaces  and  implemented  it  into  the  finite  element  program 
SOILSTRUCT-ALPHA.  The  model  is  based  on  the  Clough  and  Duncan 
(1971)  hyperbolic  formulation  and  incorporates  new  features  to  model  a 
variety  of  stress  paths.  The  Gomez-Filz-Ebeling  model  does  not  require 
any  interface  parameter  values  in  addition  to  those  used  in  the  Clough  and 
Duncan  (1971)  hyperbolic  model.  Therefore,  hyperbolic  parameter  values 
available  in  the  literature  (Clough  and  Duncan  1969,  Peterson  et  al.  1976, 
Lee  et  al.  1989)  for  a  variety  of  interfaces  can  be  used  with  the  Gomez-Filz- 
Ebeling  model. 

Gomez,  Filz,  and  Ebeling  (2000)  performed  a  series  of  interface  tests 
between  uniform,  fine  sands  and  concrete.  Some  of  their  tests  followed 
complicated  stress  paths  that  included  unloading-reloading  and  simultane¬ 
ous  changes  in  normal  and  shear  stresses.  Comparisons  between  test  results 
and  model  predictions  showed  that  the  extended  hyperbolic  model  can  pre¬ 
dict  accurately  the  response  of  interfaces  of  the  types  tested. 

Gomez,  Filz,  and  Ebeling  (2000)  also  carried  out  a  pilot-scale  lock  wall 
simulation  that  modeled  placement  and  compaction  of  backfill,  surcharge 
application,  and  changes  in  the  elevation  of  the  water  table  behind  the  wall. 
The  backfill  consisted  of  a  uniform,  fine  sand  identical  to  one  of  the  soils 
used  for  interface  testing.  The  lock  wall  simulation  test  was  also  modeled 
using  SOILSTRUCT-ALPHA.  Comparison  between  the  data  from  the  lock 
wall  simulation  and  the  results  of  the  SOILSTRUCT-ALPHA  analyses  con¬ 
firmed  the  accuracy  of  the  extended  hyperbolic  model  and  its  applicability 
for  SSI  analyses  of  lock  walls. 

The  Gomez-Filz-Ebeling  model  is  also  presumed  to  be  applicable  for 
SSI  analyses  of  multi-anchored  systems.  As  illustrated  by  the  simplified 
interface  loading  mechanism  in  Figure  1-4,  the  soil-to-wall  interface  in  a 
reinforced  concrete  slurry  wall  may  be  subject  to  simultaneous  changes  in 
shear  and  normal  stresses  as  well  as  unloading-reloading.  The  Gomez-Filz- 
Ebeling  model  is  accurate  for  predicting  the  response  of  interfaces  to  this 
type  of  loading.  While  the  model  has  been  verified  against  the  results  of 
tests  performed  on  interfaces  between  fine  sands  and  concrete,  it  has  not 
been  evaluated  against  results  of  tests  performed  on  interfaces  between 
coarse  sands  and  concrete.  In  addition,  further  work  on  interface  testing  is 
required  to  expand  the  existing  database  on  interface  hyperbolic  parameter 
values. 


1 .5  Scope  of  the  Investigation 

The  purpose  of  this  investigation  is  the  validation  of  the  Gomez-Filz- 
Ebeling  model  against  results  of  shear  tests  performed  at  the  interface 
between  a  coarse  sand  and  concrete.  Hyperbolic  parameter  values  of  the 
interface  between  the  coarse  sand  and  concrete  were  determined  from  the 
results  of  initial  loading  tests  under  constant  normal  stress.  In  addition,  an 
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interface  test  was  performed  following  a  complex  stress  path  that  included 
simultaneous  changes  in  shear  and  normal  stresses  as  well  as  unloading¬ 
reloading  of  the  interface.  The  Gomez-Filz-Ebeling  model,  with  hyperbolic 
parameter  values  determined  from  the  initial  loading  tests,  was  used  to 
estimate  the  response  of  the  interface  under  this  complex  stress  path.  The 
calculated  response  was  compared  with  the  test  measurements  to  evaluate 
the  accuracy  of  the  model. 

It  was  found  that  the  Gomez-Filz-Ebeling  model  provides  accurate  esti¬ 
mates  for  the  coarse  sand-to-concrete  interface.  It  can  be  concluded  from 
this  result  and  previous  work  that  the  Gomez-Filz-Ebeling  model  provides 
accurate  predictions  of  the  response  of  interfaces  between  concrete  and  a 
variety  of  granular  soils.  The  hyperbolic  parameters  for  the  coarse  sand-to- 
concrete  interface,  determined  from  the  tests,  also  add  to  the  database  of 
interface  parameter  values. 

The  results  of  the  laboratory  tests  performed  for  this  investigation  are 
presented  in  Chapter  2.  Results  include  those  tests  performed  to  charac¬ 
terize  the  coarse  sand  used  for  the  interface  tests,  as  well  as  from  the  inter¬ 
face  tests. 

Chapter  3  contains  a  brief  description  of  the  Gomez-Filz-Ebeling  inter¬ 
face  model.  The  hyperbolic  parameter  values  for  the  coarse  sand-to- 
concrete  interface  are  determined  and  added  to  the  database  of  interface 
parameter  values.  The  model  performance  is  evaluated  against  the  results 
of  the  interface  tests  described  in  Chapter  2. 

Chapter  4  contains  a  summary  of  the  work  performed  and  relevant  con¬ 
clusions  obtained  from  this  investigation. 

This  report  also  includes  three  appendixes.  Appendix  A  contains  the 
results  of  the  triaxial  and  consolidation  tests  performed  on  the  coarse  sand 
used  for  interface  testing.  It  also  contains  a  description  of  the  procedure  fol¬ 
lowed  for  determining  the  hyperbolic  parameter  values  of  the  coarse  sand, 
as  well  as  a  compilation  of  hyperbolic  parameter  values  for  different  soils. 
Appendix  B  contains  the  results  of  the  interface  tests  performed  and  a 
description  of  the  procedure  for  determining  the  hyperbolic  parameter  val¬ 
ues  of  the  coarse  sand-to-concrete  interface.  Finally,  all  the  symbols  and 
abbreviations  used  in  this  report  are  defined  in  Appendix  C. 
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2  Laboratory  Testing 


A  series  of  shear  tests  were  performed  on  the  interface  between  a 
coarse,  well-graded  sand  and  concrete  using  the  Large  Displacement  Shear 
Box  (LDSB).  The  objective  of  the  tests  was  to  compile  additional  experi¬ 
mental  data  on  the  response  of  soil-to-concrete  interfaces,  and  to  further 
validate  the  Gomez-Filz-Ebeling  extended  hyperbolic  model. 

This  chapter  contains  a  description  of  the  properties  of  the  sand  used  for 
this  investigation  and  results  of  the  interface  tests  performed.  Details  of  the 
techniques  for  preparation  of  the  concrete  and  soil  specimens,  and  of  the 
interface  testing  procedures,  are  given  by  Gomez,  Filz,  and  Ebeling  (2000). 


2.1  Properties  of  Blacksburg  Sand 

A  processed,  well-graded  sand,  commercially  available  for  the  prepara¬ 
tion  of  concrete  in  Blacksburg,  VA,  was  used  for  the  interface  tests.  This 
sand  is  referred  to  as  “Blacksburg  Sand”  throughout  this  report.  Table  2-1 
summarizes  the  results  of  a  series  of  laboratory  tests  performed  to  deter¬ 
mine  the  gradation,  maximum/minimum  density,  and  specific  gravity  of 
Blacksburg  Sand.  Its  grain  size  distribution  is  presented  in  Figure  2-1, 
which  shows  that  the  particle  sizes  range  from  0.2  to  5  mm.  According  to 
ASTM  D2487,  a  sand  must  have  a  value  of  the  coefficient  of  uniformity, 

C  ,  greater  than  or  equal  to  6  to  be  well  graded.  The  value  of  Cu  for  the 
Blacksburg  Sand  is  3;  therefore,  according  to  ASTM  D2487,  it  is  poorly 
graded.  However,  Blacksburg  Sand  includes  a  wider  range  of  grain  sizes 
than  did  the  sands  used  in  the  previous  study  (Gomez,  Filz,  and  Ebeling 
2000). 
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Table  2-1 

Properties  of  Blacksburg  Sand 

Parameter1 

Value 

Relevant  Standard 

Dio 

®30 

D60 

k 

0.27  mm 

0.45  mm 

0.82  mm 

3 

0.9 

ASTM  D2487 

18  kN/m3 

ASTM  D4253 

14.5  kN/m3 

ASTM  D4254 

2.65 

ASTM  D854 

1  Parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 

Examination  under  an  optical  microscope  revealed  that  the  sand  grains 
smaller  than  1  mm  are  predominantly  subangular,  with  length-to-width 
ratios  ranging  from  1  to  3.  The  larger  grains  tend  to  be  more  angular  and 
flat  (ASTM  D2488). 
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2.1.1  Triaxial  testing 


Drained  triaxial  (CD)  tests  were  performed  to  determine  the  internal  fric¬ 
tion  angle  and  hyperbolic  parameter  values  of  Blacksburg  Sand  in  a  dense 
condition.  The  specimens  for  the  tests  were  prepared  by  compaction  to  a 
relative  density  of  80  percent,  and  subjected  to  an  internal  manometric  pres¬ 
sure  of  -15  to  -20  kPa,  which  was  gradually  removed  during  application  of 
the  cell  pressure.  The  samples  were  de-aired  using  carbon  dioxide,  inun¬ 
dated  with  de-aired  distilled  water,  and  back-pressure  saturated.  The  sam¬ 
ples  in  each  set  were  consolidated  under  effective  confining  pressures 
ranging  from  69  to  276  kPa.  Shearing  was  performed  at  a  strain  rate  of 
0.25  percent  per  minute,  which  was  found  to  be  appropriate  for  pore  pres¬ 
sure  dissipation  during  previous  trials. 

The  results  of  the  tests  are  presented  graphically  in  Figure  A1  of 
Appendix  A  and  are  summarized  in  Table  2-2. 


Table  2-2 

Summary  of  Results  of  CD  Triaxial  Tests  on  Dense  Blacksburg 

Sand  (Dr  =  80%) 

Parameter1 

Value 

43.4  deg 

A<j> 

7.4  deg 

<lw 

30.7  deg 

Reference  (Figure  No.) 

Al 

1  Parameters  are  listed  and  defined  in  the  notation  (Appendix  C).  [ 

All  the  specimens  exhibited  dilation  during  shear  and  strain-softening 
after  mobilization  of  the  peak  strength.  The  peak  strength  values  measured 
during  the  tests  defined  curved  strength  envelopes.  The  value  of  the  peak 
secant  friction  angle  <|)  for  a  given  effective  confining  stress  o3'  can  be  cal¬ 
culated  from  the  following  expression  (Duncan  et  al.  1980): 


<t>=<t>0~A<t> -logic) 


(  '  \ 
£3. 

\Pa  ) 


where 


(2-1) 


<j)  =  peak  secant  friction  angle  at  a  confining  pressure  of 
101.3  kPa  (1  atm) 

A<)>  =  reduction  in  peak  secant  friction  angle  value  for  a  10-fold 
increase  in  o3' 

03'  =  minor  principal  effective  stress 
pa=  atmospheric  pressure  (101.3  kPa) 
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The  values  of  4>0  and  A<j>  were  determined  as  shown  in  Appendix  A  and 
are  given  in  Table  2-2.  The  strength  envelopes  corresponding  to  the 
strength  measured  at  a  15-percent  strain  are  linear,  and  the  corresponding 
friction  angle  <t>cv  value  is  also  given  in  Table  2-2. 


2.1.2  Consolidation  testing 

Consolidation  tests  were  performed  on  specimens  of  Blacksburg  Sand  to 
provide  additional  data  on  its  mechanical  properties.  The  specimens  were 
prepared  at  a  relative  density  of  75  percent  in  a  dry  condition.  Each  speci¬ 
men  was  consolidated  under  a  series  of  vertical  stress  increments.  The  tests 
included  a  stage  of  inundation  of  the  specimens  to  determine  the  suscepti¬ 
bility  of  Blacksburg  Sand  to  hydrocompression.  The  results  of  these  tests 
are  presented  in  Figure  A2  of  Appendix  A. 


2.1.3  Hyperbolic  parameters 

The  procedure  for  determining  the  hyperbolic  parameter  values  listed  in 
Table  2-3  is  described  in  detail  in  Appendix  A.  The  triaxial  test  data  were 
used  to  obtain  the  parameter  values  for  initial  loading,  according  to  the  pro¬ 
cedures  described  by  Duncan  et  al.  (1980).  Parameter  values  obtained  from 
the  results  from  drained  CD  triaxial  tests  on  a  variety  of  soils  are  reported 
in  Table  5  of  Duncan  et  al.  (1980).  These  values  are  reproduced  in 
Table  A2  of  Appendix  A. 


Table  2-3 

Hyperbolic  Parameter  Values  of  Dense  Blacksburg  Sand 

Hyperbolic  Parameters1 

Dense  Blacksburg  Sand 
(Dr  =  80%) 

Reference  Values  for  Soils  of 
Similar  Gradation  (from 

Duncan  et  al.  1980) 

K 

642 

340-650 

n 

0.25 

0.38-0.45 

p, 

0.55 

0.7-0.77 

Kh 

567 

230-280 

m 

0.01 

0.05-0.06 

43.4 

48-49 

A<j> 

7.4 

10-12 

1  Hyperbolic  parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 

Table  2-3  shows  ranges  of  values  for  the  modulus  exponent  n,  failure 
ratio  /?yv  bulk  modulus  exponent  m,  and  the  parameter  A<|),  which  were 
reported  by  Duncan  et  al.  (1980)  for  soils  of  similar  gradation  (Table  A2) 
Some  of  the  hyperbolic  parameter  values  for  the  dense  Blacksburg  Sand 
are  outside  the  range  of  values  obtained  from  Duncan  et  al.  (1980),  but 
they  do  not  differ  by  much. 


Chapter  2  Laboratory  Testing 


A  comparison  between  the  hyperbolic  stress-strain  relationships  calcu¬ 
lated  using  the  parameters  in  Table  2-3  and  the  data  from  the  triaxial  tests 
is  presented  in  Figure  A7  of  Appendix  A.  The  hyperbolic  model  provides  a 
good  fit  of  the  laboratory  data,  but  it  does  not  model  the  post-peak  strain¬ 
softening  behavior  or  shear-induced  dilation  of  the  soil. 


2.2  Interface  Testing 

A  set  of  tests  were  performed  at  the  interface  between  a  concrete  speci¬ 
men  and  dense  Blacksburg  Sand  using  the  LDSB.  The  purposes  of  the  tests 
were: 


a.  To  collect  additional  experimental  data  for  assessing  the  accuracy  of 
the  extended  hyperbolic  model  developed  by  Gomez,  Filz,  and 
Ebeling  (2000).  This  model  was  developed  based  on  the  results  of 
tests  performed  at  the  interface  between  concrete  and  uniform,  fine 
sands.  An  important  objective  of  the  current  research  was  to  vali¬ 
date  the  model  against  the  results  of  interface  testing  between  con¬ 
crete  and  a  well-graded,  coarse  sand,  such  as  Blacksburg  Sand. 

b.  To  add  to  the  database  on  interface  properties  available  in  the 
literature. 

For  each  of  the  tests,  the  soil  specimen  was  densified  by  vibration  to  a 
relative  density  of  80  percent,  following  the  procedures  described  by 
Gomez,  Filz,  and  Ebeling  (2000).  All  the  tests  were  performed  at  a  dis¬ 
placement  rate  of  1  mm/min  (0.04  in./min),  under  normal  pressures  ranging 
from  38  to  292  kPa. 

A  set  of  preliminary  tests  were  performed  to  study  the  influence  of  inun¬ 
dation  of  the  interface  on  the  test  results.  The  results  showed  that  inunda¬ 
tion  does  not  induce  any  significant  effect  on  the  response  of  the  interface. 
All  subsequent  tests  were  performed  in  a  dry  condition. 

Additional  details  pertaining  to  the  equipment,  specimen  preparation, 
and  testing  procedures  are  given  by  Gomez,  Filz,  and  Ebeling  (2000). 


2.2.1  Interface  testing  program 

Two  types  of  interface  tests  were  performed  for  this  report.  Four  initial 
loading  tests  were  performed  to  determine  the  peak  and  residual  friction 
angle  values  and  hyperbolic  parameter  values  of  the  interface.  These  tests 
consisted  of  shearing  the  newly  prepared  interface  under  constant  normal 
stress.  In  some  of  the  tests,  shearing  continued  until  the  residual  strength 
was  mobilized. 
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A  multidirectional  stress  path  test  was  performed  in  which  the  interface 
was  sheared  following  a  stress  path  that  involved  simultaneous  changes  in 
normal  and  shear  stresses,  as  well  as  unloading-reloading.  It  was  intended 
to  model  loading  conditions  similar  to  those  expected  at  soil-to-wall  inter¬ 
faces,  and  to  provide  a  basis  for  evaluating  the  Gomez-Filz-Ebeling  model. 
Table  2-4  summarizes  the  types  of  tests  performed  during  this  investigation 


Table  2-4 

Summary  of  Tests  Performed  at  the  Interface  Between  Concrete 
and  Dense  Blacksburg  Sand 

Normal  Stress 

Type  of  Test 

Specimen 

Test  Number 

(kPa) 

Figure  No.1 

Initial  loading 

S2001 

T2001  5 

38 

B1.B2 

(virgin  shear) 

S2002 

T2002  10 

69 

S2003 

T2003  20 

141 

S2004 

T2004_40 

292 

Multidirectional 

S2005 

T2005 

80  to  280 

B3 

stress  path 

1  Fiaures  are  included  in  Appendix  B. 

2.2.2  Results  of  interface  tests 

The  results  of  the  interface  tests  performed  for  this  investigation  are 
presented  in  the  figures  of  Appendix  B.  Table  2-4  lists  the  figures  in 
Appendix  B  that  are  relevant  to  each  interface  test. 

The  results  of  the  initial  loading  tests  are  presented  in  Figures  B1  and 
B2,  and  are  summarized  in  Table  2-5.  It  can  be  seen  that  the  peak  shear 
strength  was  mobilized  at  small  displacements  ranging  from  1  to  5  mm. 
After  mobilization  of  the  initial  peak  strength,  shear  stresses  remained 
practically  constant  before  displacement  softening  took  place. 


Table  2-5 

Summary  of  Results  of  Initial  Loading  Tests  on  the  Interface 
Between  Concrete  and  Dense  Blacksburg  Sand 

Displacement1  (mm) 

Interface  Friction  Angle1 

Interface/Soil  Friction  Angle 
Ratio 

to  Peak  A„_ 

ap 

to  Residual  As, 

Peak  5 

Residual  5r 

8/4>0 

1-5 

6-10 

32  deg 

29.5  deg 

0.74 

0.96 

1  Parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 

This  plateau  of  relatively  constant  shear  stress  may  exist  up  to  displace¬ 
ments  as  large  as  10  mm,  as  evidenced  by  test  T2003_20,  after  which  the 
residual  condition  is  attained.  In  some  cases,  a  secondary  peak  is  observed 
in  the  shear  stress-displacement  curve.  For  the  purposes  of  this  investiga¬ 
tion,  the  peak  strength  was  assumed  to  correspond  to  the  first  peak  that 
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occurred  during  the  test.  These  observations  are  consistent  with  the  results 
of  the  previous  interface  tests  reported  by  Gomez,  Filz,  and  Ebeling  (2000). 

The  peak  and  residual  interface  strength  envelopes,  which  are  shown  in 
Figure  B2  (Appendix  B),  are  both  represented  as  straight  lines  out  of  the 
origin.  As  shown  in  Figure  A1  for  the  Blacksburg  Sand  itself,  the  peak  soil 
strength  envelope  is  curved.  This  difference  between  the  curvature  of  the 
peak  soil  strength  envelope  and  the  linearity  of  the  peak  interface  strength 
envelope  was  also  observed  for  other  sands  tested  in  the  previous  study 
(Gomez,  Filz,  and  Ebeling  2000).  Figure  A1  also  shows  that  the  soil 
strength  envelope  for  the  Blacksburg  Sand  at  15  percent  strain  is  a  straight 
line  out  of  the  origin,  which  is  a  characteristic  similar  to  the  linearity  of  the 
residual  strength  of  the  interface. 

The  values  of  peak  and  residual  friction  angle,  8  and  8r,  for  the  interface 
between  Blacksburg  Sand  and  concrete  are  presented  in  Table  2-5.  These 
values  are  compared  to  the  peak  secant  friction  at  1  atmosphere  confining 
pressure,  <t>0,  and  the  15-percent  strain  friction  angle,  <|)cv,  of  the  soil.  It  can 
be  seen  that  the  values  of  8  and  8r  correspond  to  approximately  74  and 
96  percent  of  (j)0  and  (|>cv,  respectively. 

The  results  of  the  initial  loading  tests  were  used  to  determine  the  hyper¬ 
bolic  parameter  values  of  the  interface  between  concrete  and  dense 
Blacksburg  Sand,  following  the  procedure  described  in  Appendix  B. 

Figure  B3  shows  the  results  of  multidirectional  stress  path  test  T2005. 
The  test  consisted  of  the  application  of  a  complex  stress  path  that  included 
simultaneous  changes  in  shear  and  normal  stresses  and  unloading¬ 
reloading.  The  purpose  of  this  test  was  to  provide  a  basis  for  a  performance 
evaluation  of  the  extended  hyperbolic  model. 


2.3  Summary 

The  following  laboratory  activities  were  performed  for  this  investigation: 

a.  Selection  of  a  coarse  sand  specimen  for  interface  testing. 

b.  Grain  size  distribution,  minimum/maximum  density,  specific  gravity, 
consolidation  testing,  and  triaxial  testing  on  Blacksburg  Sand. 

c.  Initial  loading  tests  to  determine  the  hyperbolic  parameter  values  of 
the  interface  between  concrete  and  dense  Blacksburg  Sand. 

d.  A  multidirectional  stress  path  test  to  obtain  data  for  the  validation  of 
the  extended  hyperbolic  model. 
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A  coarse,  well-graded  sand  (Blacksburg  Sand)  was  selected  for  interface 
testing.  The  basic  properties  of  Blacksburg  Sand  were  determined  from  a 
series  of  laboratory  tests,  such  as  minimum/maximum  density  and  grain 
size  analyses.  CD  triaxial  tests  were  performed  to  determine  hyperbolic 
parameter  values  for  this  soil  (Appendix  A).  These  hyperbolic  parameter 
values  are  consistent  with  values  reported  by  Duncan  et  al.  (1980)  for  soils 
of  similar  gradation. 

Initial  loading  tests  were  performed  at  the  interface  between  concrete 
and  dense  Blacksburg  Sand  using  the  LDSB.  The  interface  hyperbolic 
parameter  values  were  determined  from  the  results  of  these  tests  as  indi¬ 
cated  in  Appendix  B. 

A  multidirectional  stress  path  test  was  performed  on  the  dense  Blacks¬ 
burg  Sand  against  concrete  interface.  The  purpose  of  this  test  was  to  obtain 
data  for  evaluating  the  accuracy  of  the  extended  hyperbolic  model  for  pre¬ 
dicting  interface  response  under  complex  stress  paths,  as  discussed  in 
Chapter  3. 
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Evaluation  of  the 
Gomez-Filz-Ebeling  Model 
for  Interfaces 


A  new  model  for  interfaces  was  developed  in  a  previous  investigation 
by  Gomez,  Filz,  and  Ebeling  (2000).  The  model  is  based  on  the  Clough  and 
Duncan  (1971)  hyperbolic  formulation,  which  was  extended  to  model  the 
response  of  interfaces  under  general  field  loading  conditions  in  lock  walls. 
The  model  captures  important  aspects  of  interface  response  while  retaining 
much  of  the  simplicity  of  the  original  Clough  and  Duncan  (1971) 
formulation. 

The  model  incorporates  three  new  elements  into  the  hyperbolic 
formulation:  (1)  a  yield  surface,  (2)  a  formulation  for  virgin  shear  under 
continuously  changing  normal  stress,  and  (3)  a  formulation  for  unloading¬ 
reloading  of  the  interface.  Gomez,  Filz,  and  Ebeling  (2000)  evaluated  the 
performance  of  the  extended  hyperbolic  model  against  the  results  of  tests 
performed  at  the  interfaces  between  concrete  and  two  uniform,  fine  sands. 
They  found  that  the  model  provides  an  accurate  estimate  of  interface 
response  for  stress  paths  such  as  those  expected  to  occur  at  the  interface 
between  a  lock  wall  and  the  backfill. 

In  this  investigation,  the  performance  of  the  extended  hyperbolic  model 
is  evaluated  against  the  results  of  tests  performed  at  the  interface  between 
coarse,  well-graded  sand  (Blacksburg  Sand)  and  concrete.  These  interface 
tests  were  described  in  Chapter  2.  It  was  found  that  predictions  of  the 
extended  hyperbolic  model  are  also  accurate  for  this  type  of  interface.  This 
suggests  that  the  model  is  applicable  for  interfaces  between  concrete  and  a 
variety  of  granular  soils.  Consequently,  it  is  possible  that  the  model  can  be 
used  for  routine  SSI  analyses  of  multi-anchored  systems  in  granular  soils, 
and  lock  walls  with  granular  backfills. 

This  chapter  contains  a  brief  review  of  the  formulation  of  the  extended 
hyperbolic  model  for  interfaces.  The  model  is  used  to  obtain  an  estimated 
response  of  the  interface  between  concrete  and  dense  Blacksburg  Sand  for 
the  complex  stress  path  applied  during  multidirectional  stress  path 
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test  T2005,  described  in  Chapter  2.  The  estimated  response  is  compared  to 
the  response  measured  during  the  test.  Additionally,  a  summary  is  pre¬ 
sented  of  hyperbolic  parameter  values  of  different  interfaces,  including 
those  reported  by  Gomez,  Filz,  and  Ebeling  (2000). 


3.1  The  Gomez-Filz-Ebeling  Model  for  Interfaces 

This  section  contains  a  summary  of  the  extended  hyperbolic  model  for 
interfaces  developed  by  Gomez,  Filz,  and  Ebeling  (2000),  who  present  a 
more  complete  description  of  the  model  formulation  and  its 
implementation. 

3.1.1  Yield  surfaces  and  loading  regions 

Based  on  analyses  of  the  results  of  staged  interface  tests,  Gomez,  Filz, 
and  Ebeling  (2000)  developed  a  formulation  for  yield  surfaces  in  interface 
shear.  According  to  their  observations,  two  yield  surfaces  can  be  defined  as 
illustrated  in  Figure  3-1.  An  upper  yield  surface  SL+  is  contained  in  the 


Figure  3-1 .  Loading  regions  for  Versions  I  and  il  of  the  Gomez-Filz-Ebeling  model 
for  interfaces 
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first  quadrant  of  the  Gn-x  plane  and  is  defined  by  the  past  maximum  shear 
stress  level  reached  during  shear.  A  lower  yield  surface  SL-  is  contained  in 
the  fourth  quadrant  and  is  defined  by  the  past  minimum  shear  stress  level. 
The  shear  stress  level  is  defined  as  the  ratio  of  shear  stress  to  shear 
strength  of  the  interface  (Clough  and  Duncan  1971).  For  an  interface  with 
zero  adhesion  intercept,  the  yield  surfaces  thus  defined  are  straight  lines 
passing  through  the  origin  of  the  Gn-x  plane. 

Two  types  of  loading  are  defined  as  illustrated  in  Figure  3-1.  Unloading¬ 
reloading  takes  place  along  any  stress  path  that  is  contained  within  the 
region  defined  by  the  SL+  and  SL-  surfaces.  Yielding  occurs  when  a  stress 
path  reaches  one  of  the  two  yield  surfaces.  As  will  be  discussed  in  Section 
3.3,  this  definition  of  loading  regions  is  used  in  Versions  I  and  II  of  the 
Gomez-Filz-Ebeling  model. 

In  some  cases  where  very  accurate  predictions  of  the  interface  response 
are  required,  it  is  convenient  to  introduce  additional  loading  regions.  Two 
transition  surfaces,  x+  and  T-,  are  illustrated  in  Figure  3-2.  They  are 
defined  by  the  past  maximum  and  past  minimum  values  of  shear  stress,  and 
consist  of  straight  lines  parallel  to  the  Gn  axis  in  the  Gn- T  plane.  The  transi¬ 
tion  regions  are  defined  as  the  areas  delimited  by  the  SL+  and  x+  surfaces 
and  by  the  SL-  and  x-  surfaces.  A  stress  path  corresponds  to  transition  load¬ 
ing  if  it  lies  within  any  of  the  two  transition  regions.  The  definition  of 


Figure  3-2.  Loading  regions  for  Version  III  of  the  G6mez-Filz-Ebeling  model  for 
interfaces 
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transition  loading  is  used  in  Version  III  of  the  G6mez-Filz-Ebeling  model 
as  discussed  in  Section  3.3.  Table  3-1  summarizes  the  definition  of  loading 
regions. 


Table  3-1 

Definition  of  Types  of  Loading  In  Interface  Shear 
(Gomez,  Filz,  and  Ebeling  2000) 

Type  of  loading 

Condition 

Yield-inducing  shear 

SL  =  SL+ 
or 

SL  =  SL-  \ 

Unloading-reloading 

SL  <  SL+  and  i  <  x+ 
and 

SL  >  SL-  and  x>x- 

Transition  loading 

SL  <  SL+  and  x>x+ 
or 

SL  >  SL-  and  r  <  x- 

3.1.2  Formulation  of  extended  hyperbolic  model  for 
interfaces  at  yield 

The  hyperbolic  model  for  interface  shear  by  Clough  and  Duncan  (1971) 
was  developed  for  initial  loading  into  the  first  quadrant  under  constant  nor¬ 
mal  stress.  Gomez,  Filz,  and  Ebeling  (2000)  extended  the  Clough  and  Dun¬ 
can  (1971)  constitutive  equation  to  loading  into  the  fourth  quadrant  as 
follows: 


X  =■ 


1 


K,-yw 


__ VA* 

( a  Y'  q  a„-tan5 


\Pa 


where 


x  =  interface  shear  stress 

As  =  displacement  along  the  interface 

Kj  =  stiffness  number 

yw  =  unit  weight  of  water 

tij  =  stiffness  exponent 

q  =  shear  direction  parameter 

R =  failure  ratio 
fl 


(3-1) 


Chapter  3  Evaluation  of  the  G6mez-Filz-Ebeling  Model  for  Interfaces 


The  shear  direction  parameter,  q,  is  equal  to  +1  if  shearing  takes  place 
into  the  first  quadrant  (x  >  0),  and  equal  to  -1  if  shearing  takes  place  into 
the  fourth  quadrant  (x  <  0). 

Equation  3-1  defines  a  three-dimensional  surface  in  the  A^-a^-x  space. 
Gomez,  Filz,  and  Ebeling  (2000)  hypothesized  that  the  interface  response 
during  shearing  along  any  yield-inducing  stress  path  is  uniquely  defined  by 
this  initial  loading  surface.  For  initial  loading  under  constant  normal 
stress,  the  following  expression  for  the  interface  shear  stiffness  Kst  is 
obtained  by  differentiation  of  Equation  3-1  with  respect  to  the  interface 
displacement,  As: 


Kst=Kr  yv 


( G  ^ 
w  n 

\Pa  J 


(l-q- Rjj  ■  SL ) 


(3-2) 


For  initial  loading  where  simultaneous  changes  in  shear  and  normal 
stresses  take  place,  i.e.  shearing  along  inclined  stress  paths,  the  interface 
stiffness,  K'sf,  is  determined  from  the  following  equation: 


Kt=Kstl 


(3-3) 


where  I  is  the  correction  factor  for  the  inclination  of  the  stress  path. 
Gomez,  Filz,  and  Ebeling  (2000)  derived  mathematically  the  following 
equation  for  the  correction  factor: 


_ l _  (3-4) 

l  +  SL-lnj^  Rfi  SL-l)-  Q’Rfj  •5'L]tan(5)  tan(0) 


where 


SL  =  stress  level 

tan(0)  =  inclination  of  the  stress  path  in  the  0n-x  plane  with  respect 
to  the  x-axis 

The  formulation  for  interface  yielding  of  the  extended  hyperbolic  model 
is  given  by  Equations  3-2,  3-3,  and  3-4.  Table  3-2  summarizes  the  state 
variables,  interface  material  properties,  and  equations  necessary  for  imple¬ 
menting  this  formulation  into  SSI  analyses.  This  formulation  does  not 
require  any  interface  properties  in  addition  to  those  defined  in  the  Clough 
and  Duncan  (1971)  hyperbolic  model. 
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Table  3-2 

Summary  of  the  Formulation  for  Interfaces  at  Yield  in  the 
Gomez-Filz-Ebeling  Model 

Conditions  for 
application 

SL  =  SL- 
or 

SL  =  SL+ 

State  variables 

SL  =  stress  level 

SL+  =  position  of  yield  surface  in  first  quadrant 

SL-  =  position  of  yield  surface  in  fourth  quadrant 
on  =  normal  stress  on  the  interface 

6  =  angle  between  the  stress  path  direction  and  the  t-axis  in  the  an-t  plane 
q=  shear  direction  parameter 

Interface 

properties 

Ki  =  stiffness  number 
rij  =  stiffness  exponent 
=  failure  ratio 
§  =  interface  friction  angle 

Equations 

k;.=i ks, 

(  \nj 

Kst  =  K,  ^  ’(I -q-Rjg^L)2 

\  Pa  ) 

/=  1 

l  +  SL-^ij(c[  RjySL-i)  -q  Rfi  tan(6) 

3.1.3  Formulation  for  unloading-reloading  of  the 
Gomez-Filz-Ebeling  model 

Three  different  versions  of  the  unload-reload  formulation  were  included 
in  the  Gomez-Filz-Ebeling  model: 

a.  Version  I.  The  interface  response  is  assumed  to  be  linear  and  normal 
stress-dependent.  No  differentiation  is  made  between  unloading¬ 
reloading  and  transition  loading.  This  version  was  intended  for  use 
in  cases  where  unloading-reloading  takes  place  along  simple  stress 
paths,  or  where  a  high  degree  of  accuracy  is  not  required  for  the 
analyses. 

b.  Version  II.  The  interface  response  is  assumed  to  be  nonlinear  and 
normal  stress-dependent.  It  follows  a  hyperbolic  shear  stress- 
displacement  relationship.  No  differentiation  is  made  between 
unloading-reloading  and  transition  loading.  This  formulation  is 
simple  to  implement  in  finite  element  programs,  yet  it  is  accurate 
for  most  loading  conditions. 
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c.  Version  III.  As  in  Version  II,  the  interface  response  is  assumed  to  be 
nonlinear  and  normal  stress-dependent.  However,  the  mathematical 
expressions  used  to  convey  this  nonlinear  response  are  more  com¬ 
plex  than  those  of  Version  II.  Version  III  is  the  most  accurate  ver¬ 
sion  of  the  model.  It  is  also  the  most  complex  because  it  introduces 
additional  sets  of  state  variables  that  are  needed  to  differentiate 
between  unloading-reloading  and  transition  loading. 

In  all  three  versions  of  the  unload-reload  formulation,  the  interface 
response  is  independent  of  the  orientation  of  the  stress  path.  In  SSI  analy¬ 
ses,  any  of  the  versions  of  the  unload-reload  formulation  can  be  applied  in 
conjunction  with  the  formulation  for  interfaces  at  yield  introduced  pre¬ 
viously.  The  following  sections  describe  the  mathematical  formulation  of 
the  three  versions  of  the  unload-reload  formulation. 

3.1.3.1  Version  I.  In  this  version,  the  interface  response  is  assumed  lin¬ 
ear  and  normal  stress-dependent  in  both  unloading-reloading  and  transition 
regions.  The  interface  tangent  stiffness  is  determined  from  the  following 
expression: 


K  =  Kurj  -yw 


ra  YV 

n 


Pa) 


(3-5) 


where  Kur ■  is  the  interface  stiffness  number  for  unloading-reloading. 

Table  3-3  contains  a  summary  of  this  version  of  the  model.  The  formula¬ 
tion  is  applicable  if  the  stress-state  is  inside  the  region  limited  by  the  yield 
surfaces  SL+  and  SL-,  in  the  on-x  plane.  There  are  no  new  state  variables  to 
be  considered  for  this  formulation  in  addition  to  those  listed  in  Table  3-2. 
The  application  of  the  model  is  very  simple  once  a  value  of  Kurj  is 
selected.  It  is  noted  that  the  interface  response  during  unloading-reloading 
is  assumed  to  be  independent  of  the  inclination  of  the  stress  path. 

This  formulation  is  appropriate  for  analyses  that  do  not  require  a  high 
degree  of  accuracy.  In  cases  where  the  interface  may  be  subjected  to  com¬ 
plex  unloading-reloading  paths,  a  single  value  of  Kurj  may  not  provide  an 
accurate  approximation  to  the  nonlinear  response  of  the  interface.  Version 
II,  described  in  the  following  section,  accounts  for  nonlinearity  of  the  inter¬ 
face  response  during  unloading-reloading. 

3.1.3.2  Version  II.  In  Version  II,  the  tangent  interface  stiffness  during 
unloading-reloading  is  determined  from  the  following  expression: 


K  =  Kn  -j 


'a' 


\Pa) 


(3-6) 
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Table  3-3 

Formulation  for  Unloading-Reloading  and  Transition  Loading  In  the  Gomez-Fllz-Ebeling  Model  for  Interfaces 
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Note:  Symbols  are  defined  in  the  notation  (Appendix  C). 


where 


SLo  =  shear  stress  level  at  the  origin 
a  =  dimensionless  scaling  factor  for  unloading-reloading 

The  origin  of  unloading-reloading  is  defined  as  the  last  point  where  a 
change  in  the  direction  of  shear  took  place.  The  stress  level  at  the  origin 
SLa  is  calculated  from  the  following  expression: 


SL0  = 


•tan(8) 


(3-7) 


where 


%o  =  shear  stress  at  the  origin 
o  =  normal  stress  at  the  origin 

no 

According  to  this  formulation,  the  interface  response  during  unloading¬ 
reloading  is  characterized  by  a  hyperbolic  shear  stress-displacement  rela¬ 
tionship  that  is  qualitatively  identical  to  that  assumed  for  yield-inducing 
shear.  However,  the  interface  response  during  unloading-reloading  is 
assumed  to  be  independent  of  the  inclination  of  the  stress  path.  In  Ver¬ 
sion  II,  Expression  3-6  is  assumed  applicable  to  both  unloading-reloading 
and  transition  loading. 

As  can  be  seen  in  Table  3-3,  this  version  of  the  model  introduces  two 
new  interface  properties:  the  stiffness  number  for  unloading-reloading, 

K  ,  and  the  scaling  factor,  a.  These  property  values  can  be  determined 
through  a  trial-and-error  process  from  the  results  of  unload-reload  interface 
tests.  They  can  also  be  related  to  other  interface  properties  without  the 
need  for  their  experimental  determination.  A  procedure  for  the  determina¬ 
tion  of  a  and  Kurj  from  other  interface  properties  was  proposed  by  Gomez, 
Filz,  and  Ebeling  (2000).  This  procedure  is  summarized  in  Section  3.1.4 
(Determination  of  the  model  parameter  values). 

3.1.3.3  Version  III.  This  version  is  intended  for  use  in  analyses  where 
very  accurate  modeling  of  the  interface  response  is  required.  It  is  a  more 
sophisticated  formulation  that  accounts  for  the  type  of  interface  response 
observed  in  the  transition  region.  From  a  practical  point  of  view,  it  is  not 
known  whether  Version  III  introduces  significant  gains  in  accuracy  for  the 
modeling  of  soil-to-wall  interfaces  with  respect  to  Versions  I  and  II.  Addi¬ 
tional  studies  would  be  necessary  to  investigate  the  comparative  advan¬ 
tages  of  Version  III  for  practical  problems. 
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In  Version  III,  the  interface  response  during  unloading-reloading  is  iden¬ 
tical  to  Version  II.  Therefore,  for  stress  paths  contained  inside  the 
unloading-reloading  region  defined  in  Figure  3-2,  the  tangent  interface  stiff¬ 
ness  can  be  determined  from  Equation  3-6. 

The  difference  between  the  formulations  of  Versions  II  and  III  corre¬ 
sponds  to  the  transition  region.  In  Version  II,  the  formulations  for  the 
unloading-reloading  and  the  transition  regions  are  identical  (Equations  3-6 
and  3-7).  In  Version  III,  the  interface  stiffness  at  any  point  inside  the  transi¬ 
tion  region  is  calculated  using  the  following  expression: 


jy-f  _  ts  tS 

A St  "  A5?2 


mk(SL-Sl!s} 

10  V  '-Yh  ■ 


where 


K  ts  =  transition  stiffness  number 
sn 

mk  =  stiffness  degradation  parameter 
SLts  =  transition  stress  level 


(3-8) 


Equation  3-8  represents  the  formulation  for  transition  loading  in  Version 
III  of  the  G6mez-Filz-Ebeling  model.  The  transition  stiffness  number  Ksnts, 
stiffness  degradation  parameter  mk,  and  transition  stress  level  SLts  are  three 
of  the  additional  state  variables  required  for  implementation  of  this 
formulation. 


The  transition  stiffness  number  can  be  determined  from  the  following 
expression: 
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where 

K'  ts  =  interface  stiffness  at  the  point  where  the  stress  path  intersects 
the  transition  surface  upon  entering  the  corresponding 
transition  region 

a  ts  =  normal  stress  at  the  point  where  the  stress  path  intersects  the 
transition  surface  upon  entering  the  corresponding  transition 
region 
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The  stiffness  degradation  parameter,  mk,  is  determined  from  the 
following  expression: 


log(^„W)-log(^„yi) 

SLts-SLys 


(3-10) 


where 

K  ys  =  yield  stiffness  number 

sn  J 

SLys  =  stress  level  for  the  current  position  of  the  yield  surface 

The  yield  stiffness  number,  Ksnys,  is  the  value  of  normalized  stiffness 
corresponding  to  the  current  position  of  the  yield  surface,  SLys,  and  is 
calculated  from  the  following  equation: 

Kmys  =  KrQ.-q-RjS-Sl7s'f  (3-11) 

The  current  position  of  the  yield  surface,  SUS,  is  defined  mathemati¬ 
cally  as 

SLys  =  SL+  for  shear  inside  the  transition  zone  of  the  first  quadrant 
=  SL-  for  shear  inside  the  transition  zone  of  the  fourth  quadrant 

(3-12) 


The  value  of  the  transition  stress  level,  SLts,  is  determined  from  the 
following  expression: 


SLts 


onts  tan(8  ) 


(3-13) 


where  Xts  is  the  shear  stress  at  the  point  where  the  stress  path  intersects 
the  transition  surface  upon  entering  the  corresponding  transition  region. 

Table  3-3  summarizes  the  formulation  of  Version  III.  This  version  of  the 
model  does  not  introduce  any  new  interface  properties  with  respect  to  those 
required  in  Version  II.  A  total  of  five  additional  state  variables  are  required 
for  Version  III. 
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3.1.4  Determination  of  the  model  parameter  values 

The  application  of  the  extended  hyperbolic  model  requires  the  determi¬ 
nation  of  the  values  of  the  interface  parameters  listed  in  Table  3-4.  Gomez, 
Filz,  and  Ebeling  (2000)  presented  a  detailed  description  of  the  procedure 
for  the  determination  of  these  parameter  values,  which  is  summarized  in 
the  following  paragraphs. 


Table  3-4 

Summary  of  Material  Properties  Required  In  the 

Gomez-Filz-Ebeling  Model  for  Interfaces 

Type  of  Loading 

Formulation 

Parameters1 

Yield-inducing  shear 

Formulation  for  yield-inducing 
shear 

Ki  =  stiffness  number 
r)j=  stiffness  exponent 

Rtj=  failure  ratio 

S  =  interface  friction  angle 

Unloading-reloading 

Version  1 

Kurj=  stiffness  number  for 
unloading-reloading 
ny=  stiffness  exponent 

8  =  interface  friction  angle 

Version  II 

KUrj=  stiffness  number  for 
unloading-reloading 

Dj=  stiffness  exponent 

8  =  interface  friction  angle 

Rij  =  failure  ratio 
a  =  scaling  factor 

Transition  loading 

Version  III 

Kurj=  stiffness  number  for 
unloading-reloading 
nj  =  stiffness  exponent 

8  =  interface  friction  angle 

Rtj  =  failure  ratio 
a  =  scaling  factor 

1  The  value  of  each  of  the  parameters  n;,  5,  and  Rfj  is  identical  for  yield-inducing  shear, 
unloading-reloading,  and  transition  loading. 

3.1.4.1  Determination  of  parameter  values  for  interfaces  at  yield. 

The  implementation  of  the  formulation  for  yield-inducing  shear  of  the 
extended  hyperbolic  model  does  not  require  any  parameters  in  addition  to 
those  considered  in  the  Clough  and  Duncan  (1971)  hyperbolic  model.  A 
detailed  description  of  the  procedure  for  determining  the  hyperbolic 
parameter  values  is  presented  in  Appendix  B. 

Gomez,  Filz,  and  Ebeling  (2000)  observed  that  the  simple  procedure  for 
determining  the  hyperbolic  parameter  values  developed  by  Clough  and 
Duncan  (1971)  is  accurate  enough  for  implementation  of  the  hyperbolic 
model  in  cases  of  shear  under  constant  normal  stress.  However,  for  shear 
along  inclined  stress  paths,  it  is  convenient  to  verify  and,  if  necessary,  to 
adjust  the  hyperbolic  parameter  values  before  application  of  the  extended 
hyperbolic  model.  The  procedure  for  adjusting  the  hyperbolic  parameters  is 
based  on  normalization  of  the  test  data.  This  procedure  is  summarized  in 
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Appendix  B.  Further  details  on  normalization  of  interface  test  data  and 
adjustment  of  the  hyperbolic  parameters  can  be  found  in  Gomez,  Filz,  and 
Ebeling  (2000). 

3. 1.4.2  Determination  of  interface  parameter  values  for  unloading¬ 
reloading.  The  stiffness  number  for  unloading-reloading,  KUrj ,  is  required 
for  Versions  I,  II,  and  III  of  the  model.  The  value  of  KUrj  can  be  determined 
experimentally  from  the  results  of  unload-reload  tests,  or  can  be  estimated 
using  the  following  equation: 

Kurj  =  CkKj  (3-14) 

where  Ck  is  the  interface  stiffness  ratio.  The  value  of  the  stiffness  ratio 
can  be  estimated  according  to  the  expressions  presented  in  Table  3-5. 


Table  3-5 

Recommended  Values  of  the  Stiffness  Ratio  Ck 
(Gomez,  Filz,  and  Ebeling  2000) 

Version  of  the  Model 

Correlation  with  Rfj 

Expected  Range 

1 

Ck  =  0.5  ■  (1  +  Rfjj2 

1  to  2 

II  and  III 

2  to  4 

The  scaling  factor,  a,  is  used  in  Versions  II  and  III  of  the  model.  Its 
value  can  also  be  determined  experimentally  from  the  results  of  unload- 
reload  tests,  or  it  can  be  estimated  according  to  the  expression  presented  in 
Table  3-6. 


Table  3-6 

Recommended  Values  of  the  Scaling  Factor  a 
(Gomez,  Filz,  and  Ebeling  2000) 

Version  of  the  Model 

Correlation  with  Rf 

Expected  Range 

1 

Not  applicable 

— 

II  and  III 

01=1  +Rn 

1 .4  to  2 

3.2  Interface  Hyperbolic  Parameter  Values 

Hyperbolic  parameter  values  for  the  interface  between  dense  Blacksburg 
Sand  and  concrete  were  determined  from  the  results  of  the  initial  loading 
tests  described  in  Chapter  2,  as  described  in  Appendix  B.  These  values  are 
summarized  in  Table  3-7. 
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Table  3-7 

Hyperbolic  Parameter  Values  for  the  Interface  Between  Dense 
Blacksburg  Sand  and  Concrete 


Parameter 

Parameter  Value 

K, 

23000 

n, 

0.8 

0.76 

8 

31 .6  deg 

The  parameter  values  for  the  interface  between  dense  Blacksburg  Sand 
and  concrete  are  reproduced  in  Table  3-8,  together  with  the  parameter  val¬ 
ues  reported  in  the  literature  for  other  interfaces.  This  table  may  be  useful 
for  the  selection  of  parameter  values  in  cases  where  results  of  initial  load¬ 
ing  tests  are  not  available. 


Table  3-8 

Summary  of  Hyperbolic  Parameter  Values  for  Interfaces  between  Concrete  and 
Different  Types  of  Soils 


Soli 

Concrete  Surface 

Relative  Density  of 
Soli 

■ 

■ 

Well-graded  Blacksburg  Sand 

Cast  against  plywood 

80%  (dense) 

23000 

0.8 

0.76 

31.6 

Uniform  density  sand  (Gomez,  Filz, 
and  Ebeling  2000) 

Cast  against  plywood 

49%  (medium  dense) 

21800 

0.81 

0.88 

29.3 

75%  (dense) 

26600 

0.83 

0.85 

31 

Uniform  Light  Castle  Sand  (Gomez, 
Filz,  and  Ebeling  2000) 

Cast  against  plywood 

80%  (dense) 

20700 

0.79 

Uniform  sand  (Peterson  et  al.  1976) 

Smooth  (cast  against 
glass  plate) 

Loose 

10200 

0.87 

0.71 

30 

Medium  dense 

12700 

0.84 

0.62 

28.7 

Dense 

8400 

1.17 

0.40 

31 .2 

Rough 

Loose 

10000 

0.83 

0.85 

31 .6 

Medium  dense 

11900 

0.71 

0.78 

29.8 

Dense 

10400 

0.70 

0.41 

32.9 

Well-graded  sand  (Peterson  et  al. 
1976) 

Smooth  (cast  against 
glass  plate) 

Loose 

12000 

Qiisi 

Medium  dense 

9200 

Dense 

10500 

Rough 

Loose 

7700 

Medium  dense 

13100 

HPIMI 

Dense 

14800 

Poorly  graded  sand  (Clough  and 

Duncan  1 969) 

Smooth 

Dense 

Ottawa  Sand  50-60  (Lee  et  al.  1989) 

Smooth 

Dense 

Rough 
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There  are  a  larger  number  of  reference  values  of  hyperbolic  parameters 
for  soils  than  there  are  for  interfaces.  If  neither  interface  test  results  nor 
reference  parameters  are  available,  it  may  be  possible  to  relate  the  interface 
parameter  values  to  the  hyperbolic  parameter  values  of  the  soil.  Table  3-9 
compares  the  hyperbolic  parameter  values  for  the  interface  between  dense 
Blacksburg  Sand  and  concrete,  and  for  the  three  interfaces  tested  by 
Gomez,  Filz,  and  Ebeling  (2000),  to  the  hyperbolic  parameter  values 
obtained  from  CD  triaxial  tests  on  isotropically  consolidated  soil 
specimens. 


Table  3-9 

Comparison  of  the  Hyperbolic  Parameter  Values  of  Several 
Interfaces  to  Hyperbolic  Parameter  Values  of  the  Soils 

Type  of  Interface 

Hyperbolic  Parameter  Values 

Interface-to-Soil 
Parameter  Ratio 

Interface 

Soil 

Dense  Blacksburg  Sand 

Ki=  23000 

K=  642 

36 

rij  =  0.8 

n=  0.25 

3.2 

Rfj  =  0.76 

Rt=  0.55 

1.37 

8  =  31 .6  deg 

<|>o  =  43.4  deg 

0.73 

Dense  density 

Ki=  26600 

/f  =  870 

31 

sand-to-concrete 

rij  =  0.83 

n  =  0.63 

1.32 

(Gomez,  Filz,  and 

Rfj  =  0.85 

Rf=  0.85 

0.99 

Ebeling  2000) 

8  =  31  deg 

(|>o  =  36  deg 

0.86 

Medium  density 

Ki=  21900 

K=  1200 

18 

sand-to-concrete 

ny  =  0.83 

n  =  0.63 

1.32 

(Gomez,  Filz,  and 

Rfj  =  0.85 

Rf=  0.85 

1.03 

Ebeling  2000) 

8  =  31  deg 

4>o  =  41  deg 

0.71 

Dense  Light  Castle 

Ki=  20700 

/C=  690 

30 

sand-to-concrete 

nj  =  0.79 

n  =  0.79 

1.0 

(Gomez,  Filz,  and 

Rfj  =  0.79 

Rt  =  0.81 

0.97 

Ebeling  2000) 

8  =  33.7  deg 

<(>o  =  42.4  deg 

0.79 

Ratios  of  interface-to-soil  hyperbolic  parameter  values  are  presented  in 
Table  3-9.  These  ratios  are  based  on  a  limited  database,  and  they  should  not 
be  used  as  a  reference  for  interfaces  that  differ  substantially  from  the  inter¬ 
faces  tested  during  this  investigation. 

Table  3-10  summarizes  the  values  of  parameters  Kurj  and  a  for  the  inter¬ 
face  between  dense  Blacksburg  Sand  and  for  other  types  of  interfaces. 
These  are  estimated  values,  determined  using  the  criteria  presented  in 
Tables  3-5  and  3-6,  the  hyperbolic  parameter  values  listed  in  Table  3-8,  and 
Equation  3-14. 
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Table  3-10 

Summary  of  Values  of  Model  Parameters  Kurj  and  a  for  Several  Interfaces 


Versions  II  and  III 


Well-graded  Blacksburg  Sand 


Concrete  Surface  Soil 


Cast  against 
plywood 


Relative  Density  of  Version 


Uniform  density  sand  (G6mez,  Filz,  Cast  against 


and  Ebeling  2000) 


plywood 


Uniform  Light  Castle  Sand  (G6mez,  Cast  against 
Filz,  and  Ebeling  2000)  plywood 


Uniform  sand  (Peterson  et  al.  1 976)  Smooth  (cast 

against  glass 


Well-graded  sand  (Peterson  et  al. 
1976) 


Smooth  (cast 
against  glass 
plate) 


Poorly  graded  sand  (Clough  and 
Duncan  1 969) 


Smooth 


Ottawa  Sand  50-60  (Lee  et  al.  1 989)  Smooth 


Medium  dense 


Medium  dense 


131100 


262260 


The  values  of  parameters  Kurj  and  a  are  used  for  modeling  unloading-reloading  and  transition  loading  in  the  Gomez-Filz-Ebeling 
model. 

The  parameter  values  presented  in  this  table  were  estimated  using  the  criteria  in  Tables  3-5  and  3-6,  the  parameter  values  for  Kt 
and  Rf,  in  Table  3-8,  and  Equation  3-14. 
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3.3  Evaluation  of  the  Extended  Hyperbolic  Model 


Gomez,  Filz,  and  Ebeling  (2000)  carried  out  comparisons  between  the 
interface  response  calculated  using  the  extended  hyperbolic  model  and  the 
results  of  interface  tests.  They  found  that  the  model  was  accurate  to  predict 
the  response  of  interfaces  between  uniform,  fine  sands  and  concrete  sub¬ 
jected  to  complex  stress  paths.  In  this  section,  further  evaluation  of  the 
model  is  carried  out  using  the  results  of  interface  test  T2005  performed  dur¬ 
ing  this  investigation  using  the  coarse,  well-graded  Blacksburg  Sand. 

Figure  3-3a  shows  the  stress  path  applied  during  interface  test  T2005. 
The  following  types  of  loading  can  be  identified  on  this  stress  path: 

a.  Initial  loading  to  point  A:  yield  inducing-shear. 

b.  Path  A-B:  yield-inducing  shear. 

c.  Path  B-C-Tl:  unloading-reloading. 

d.  Path  Tl-Yl:  transition  loading. 

e.  Path  Yl-D:  yield-inducing  shear. 

/.  Path  D-E :  unloading. 

g.  Path  E-F:  yield-inducing  shear. 

h.  Path  F-Y2:  unloading-reloading. 

i.  Path  Y2-G-H :  yield-inducing  shear  to  failure. 

The  shear  stress-displacement  response  was  estimated  using  all  versions 
of  the  Gomez-Filz-Ebeling  model.  The  interface  parameter  values  used  are 
those  listed  in  Table  3-7.  Comparisons  between  the  measured  interface 
response  and  the  results  obtained  from  Versions  I,  II,  and  III  are  presented 
in  Figures  3-3b,  3-3c,  and  3-3d,  respectively. 

All  three  versions  of  the  model  provide  good  approximations  of  the 
interface  response  for  this  stress  path.  The  most  significant  differences 
between  the  three  versions  correspond  to  the  calculated  interface  response 
during  transition  loading  from  77  to  Yl,  and  during  unloading-reloading 
from  D  to  E  and  from  F  to  Y2.  The  interface  response  from  T1  to  Yl,  esti¬ 
mated  from  Versions  I  and  II,  is  stiffer  than  the  measured  interface 
response.  A  sharp  break  in  the  shear  stress-displacement  curve  is  obtained 
at  point  Yl  from  the  model  calculations,  which  does  not  fit  the  observed 
response.  Version  III  provides  a  more  accurate  estimate  of  the  interface 
stiffness  during  transition  loading,  and  a  smoother  transition  from  transi¬ 
tion  loading  to  initial  loading  at  point  Yl. 
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stress,  x  (kPa) 
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0.00  0.25  0.50  0.75  1.00  1.25  1.50  1.75  2.00  2.25  2.50 

Interface  displacement,  As  (mm) 

c.  Comparison  between  Version  II  for  unloading-reloading  and  test  data 
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Interface  displacement,  As  (mm) 


d.  Comparison  between  Version  III  for  transition  loading  and  test  data 


Figure  3-3.  (Concluded) 
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The  interface  response  estimated  using  Version  I  for  unloading¬ 
reloading  from  D  to  E  and  from  F  to  Y2  is  linear  and  does  not  model  the 
interface  response  accurately.  Versions  II  and  III  capture  the  nonlinearity  of 
the  interface  response. 

It  can  be  seen  that  the  formulation  for  yield-inducing  shear  provides 
accurate  estimates  of  the  interface  response  in  the  first  and  fourth  quad¬ 
rants  for  all  directions  of  the  stress  path  applied  during  this  test.  It  is  impor¬ 
tant  to  note  that,  during  shearing  between  points  G  and  H,  the  measured 
interface  response  shows  stiffness  values  that  are  lower  than  zero.  That  is, 
interface  displacement  increases  as  shear  stresses  decrease.  This  negative 
stiffness  is  not  due  to  post-peak  softening  of  the  interface  because  it  takes 
place  before  the  stress  path  intersects  the  interface  strength  envelope.  Nega¬ 
tive  values  of  interface  stiffness  can  occur  for  certain  inclinations  of  the 
stress  path  (see  Figure  4-16  in  Gomez,  Filz,  and  Ebeling  2000).  As  can  be 
seen  in  Figure  3-3,  the  formulation  for  interfaces  at  yield  of  the  extended 
hyperbolic  model  captures  this  feature  of  the  interface  response. 

For  shearing  beyond  point  H  (Figure  3-3),  once  failure  of  the  interface 
is  attained,  the  model  does  not  provide  accurate  estimates  of  the  interface 
response  because  the  stiffness  is  assumed  to  be  zero  at  failure  in  the  model. 

In  summary,  comparison  of  the  model  calculations  to  the  results  of  inter¬ 
face  test  T2005  confirms  the  accuracy  of  the  extended  hyperbolic  model  in 
the  prediction  of  the  response  of  sand-to-concrete  interfaces  subjected  to  a 
variety  of  complex  stress  paths. 


3.4  Input  of  Model  Parameters  in 
SOILSTRUCT-ALPHA 


The  formulation  for  yield-inducing  shear  and  Version  II  of  the 
formulation  for  unloading-reloading  of  the  Gomez-Filz-Ebeling  model 
were  implemented  into  the  finite  element  program  SOILSTRUCT-ALPHA. 
Table  3-11  summarizes  the  soil  and  interface  parameters  required  as  input 
for  SOILSTRUCT-ALPHA  analyses. 
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Table  3-11 

Summary  of  Soils  and  Interface  Material  Properties  Required  as 
Input  for  SOILSTRUCT-ALPHA  Analyses 

Materials 

Input  Parameters1 

Soils 

K  =  modulus  number 
n  =  modulus  exponent 

Kur-  modulus  number  for  unloading-reloading 

Rt=  failure  ratio 
(ji  =  internal  friction  angle 
vnom  =  nominal  Poisson's  ratio 

K0  =  at-rest  earth  pressure  coefficient 
ymoist,  "/sat  =  moist  and  saturated  unit  weight,  respectively 

Interface 

Ki  =  stiffness  number 
rij  =  stiffness  exponent 

Rtj  =  failure  ratio 

S  =  interface  friction  angle 

1  The  values  of  interface  parameters  K  and  a  for  unloading-reloading  are  not  input  in 
SOILSTRUCT-ALPHA.  They  are  estimated  internally  using  the  criteria  in  Tables  3-5  and  3-6,  the 

input  values  for  K,  and  Rfl,  and  Equation  3-14. 

3.4.1  Soil  parameter  input 

The  values  of  parameters  K,  n,  Kur,  Rp  and  <j>  for  soils  can  be  determined 
experimentally  from  the  results  of  triaxial  tests,  or  they  can  be  estimated 
using  values  reported  in  the  literature. 

Values  of  the  unload-reload  modulus  number,  Kur,  for  soils  are  not  read¬ 
ily  available  in  the  literature.  According  to  Duncan  et  al.  (1980),  the  value 
of  Kur  is  always  larger  than  the  value  of  the  modulus  number  K  for  primary 
loading.  For  stiff  soils  such  as  dense  sands,  Kur  may  be  20  percent  greater 
than  K.  For  soft  soils  such  as  loose  sands,  Kur  may  be  three  times  as  large 
as  K.  The  value  of  the  modulus  exponent  n  is  always  very  similar  for  pri¬ 
mary  loading  and  unloading,  and  is  often  assumed  to  be  the  same. 

Many  granular  soils  present  strength  envelopes  that  are  curvilinear  and 
cannot  be  defined  by  a  single  value  of  cj>.  In  these  cases,  a  value  of  secant 
friction  angle  <|)  can  be  assigned  to  the  soil  that  corresponds  to  the  range  of 
confining  stresses  expected  to  occur.  This  value  of  secant  friction  angle  can 
be  calculated  using  the  following  equation: 


<t>=<t>o-A<Hog10 


\Pa) 


(3-1  bis) 


Table  A2  summarizes  values  of  K,  n,  Rp  <|>0  ,  and  Acj)  for  several  types  of 
soils  at  varying  relative  densities.  These  values  were  extracted  from 
Table  5  in  Duncan  et  al.  (1980).  They  can  be  used  as  a  basis  for  selection 
of  parameters  when  triaxial  test  results  are  not  available,  or  for  confirma¬ 
tion  of  the  experimental  parameter  values. 
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Values  for  the  nominal  Poisson’s  ratio  \npm  are  determined  from  a  one¬ 
dimensional  (1-D)  column  analysis  as  described  by  Ebeling  and  Wahl 
(1997).  In  the  1-D  column  analysis,  the  horizontal  earth  pressure  coeffi¬ 
cient  Kh  resulting  from  incremental  backfill  placement  of  the  soil  column 
is  related  to  vnom,  as  described  in  Ebeling  and  Wahl  (1997). 

The  at-rest  earth  pressure  coefficient,  K0,  input  as  a  material  property  in 
SOILSTRUCT- ALPHA  is  assigned  only  to  newly  placed  fill  lifts.  Because 
newly  placed  lifts  are  usually  assumed  to  behave  as  a  dense  fluid,  the  value 
of  K  for  input  in  SOILSTRUCT- ALPHA  is  often  assumed  equal  to  1. 


3.4.2  Interface  parameter  Input 


The  values  of  parameters  Kj,  rij,  R^,  and  8  for  soil-to-structure  interfaces 
can  be  determined  experimentally,  or  they  can  be  estimated  using  values 
available  in  the  literature.  Table  3-8  summarizes  parameter  values  for  sev¬ 
eral  types  of  interfaces.  These  values  can  be  used  as  the  basis  for  assigning 
interface  properties  when  test  results  are  not  available,  or  for  verification 
of  experimental  parameter  values. 

It  must  be  noted  that,  in  SOILSTRUCT-ALPHA  analyses,  it  is  not  neces¬ 
sary  to  introduce  the  values  for  unload-reload  parameters  Kurj  and  a.  They 
are  estimated  internally  using  the  criteria  in  Tables  3-5  and  3-6,  the  input 
values  for  Kj  and  R^,  and  Equation  3-14. 


3.5  Summary  and  Conclusions 

The  extended  hyperbolic  model  developed  by  G6mez,  Filz,  and  Ebeling 
(2000)  was  evaluated  against  the  results  of  shear  tests  performed  at  the 
interface  between  a  coarse,  well-graded  sand  and  concrete.  It  was  con¬ 
firmed  that  the  model  can  predict  accurately  the  response  of  interfaces 
between  sand  and  concrete.  Principal  advantages  of  the  model  are  these: 


a.  It  has  a  simple  mathematical  formulation. 

b.  Hyperbolic  parameter  values  for  different  types  of  interfaces  are 
available  in  the  literature. 

c.  It  captures  the  main  features  of  the  interface  response  under  simulta¬ 
neous  changes  in  shear  and  normal  stress  and  during  unloading¬ 
reloading. 

d.  It  provides  accurate  estimates  of  the  response  of  a  variety  of  soil-to- 
concrete  interfaces  subjected  to  complex  stress  paths. 

e.  It  is  relatively  easy  to  implement  in  SSI  analyses. 
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f.  It  establishes  a  framework  for  future  work  on  plasticity-based  inter¬ 
face  models. 

The  formulations  for  yield-inducing  shear,  and  for  unloading-reloading 
Version  II  of  the  extended  hyperbolic  model,  were  implemented  by  Gomez, 
Filz,  and  Ebeling  (2000)  in  the  finite  element  program  SOILSTRUCT- 
ALPHA.  This  program  is  commonly  used  by  the  Corps  of  Engineers  for 
analyses  of  earth-retaining  structures.  Through  finite  element  analyses  of  a 
pilot-scale  lock  wall  simulation,  Gomez,  Filz,  and  Ebeling  (2000)  found 
that  these  formulations  could  be  used  for  predicting  vertical  shear  forces  in 
retaining  structures. 

The  model  has  several  limitations: 

a.  It  does  not  model  displacements  normal  to  the  interface,  and  the  inter¬ 
face  thickness  is  implicitly  assumed  as  zero.  Consequently,  in  finite 
element  analyses,  a  large  normal  stiffness  must  be  assigned  to  inter¬ 
face  elements  in  order  to  reduce  overlapping  of  adjacent  two- 
dimensional  elements.  In  addition,  it  cannot  model  the  generation  of 
normal  stresses  due  to  restrained  dilation  of  the  interface  during 
shear  between  two  stiff,  rough  media.  This  may  not  be  important  for 
analyses  of  stiff  structures  retaining  relatively  compressible  soils. 

b.  It  does  not  model  displacement  softening  of  the  interface.  According 
to  the  experimental  data  collected  during  this  investigation,  displace¬ 
ment  softening  may  take  place  in  interfaces  subjected  to  relative 
displacements  of  6  to  20  mm.  Therefore,  in  cases  where  larger  mag¬ 
nitudes  of  interface  displacement  take  place,  the  model  cannot  pro¬ 
vide  accurate  predictions  of  the  displacement-softening  interface 
response. 

c.  The  model  predicts  interface  stiffness  values  that  are  zero  or  negative 
for  certain  loading  combinations.  For  implementation  of  the  model 

in  finite  element  programs,  it  is  then  necessary  to  use  appropriate 
stiffness  values  and  numerical  procedures  (see  Stankowsky, 

Runesson,  and  Sture  1993,  for  example)  in  order  to  prevent  numeri¬ 
cal  problems.  It  was  found  during  this  investigation  that  negative 
stiffness  values  do  occur  for  certain  stress  path  inclinations.  How¬ 
ever,  if  finite  element  analyses  of  multi-anchored  systems  show  that 
these  types  of  loading  are  common,  it  may  be  necessary  to  perform 
additional  experimental  work  to  verify  the  accuracy  of  the  model 
predictions. 

d.  The  model  has  not  yet  been  used  for  routine  analyses  of  multi- 
anchored  systems.  However,  it  was  used  successfully  by  Gomez, 

Filz,  and  Ebeling  (2000)  to  estimate  vertical  shear  forces  in  a  pilot- 
scale  lock  wall  simulation. 
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4  Summary  and  Conclusions 


Soil-to-structure  interfaces  in  multi-anchored  systems  may  be  subjected 
to  complex  stress  paths,  which  may  include  simultaneous  changes  in  shear 
and  normal  stress  and  unloading-reloading.  Accurate  modeling  of  the  inter¬ 
face  response  to  this  type  of  loading  may  be  required  to  develop  improved 
design  procedures  of  multi-anchored  systems.  Gomez,  Filz,  and  Ebeling 
(2000)  developed  the  extended  hyperbolic  model  for  interfaces.  Based  on 
comparisons  between  model  predictions  and  results  of  laboratory  interface 
tests,  they  showed  that  the  model  provides  accurate  estimates  of  the 
response  of  interfaces  between  fine  sands  and  concrete  subjected  to  a  vari¬ 
ety  of  stress  paths. 

G6mez,  Filz,  and  Ebeling  (2000)  implemented  the  extended  hyperbolic 
model  in  the  finite  element  program  SOILSTRUCT-ALPHA,  which  is 
commonly  used  for  analyses  of  Corps  of  Engineers  structures,  such  as  lock 
walls.  They  also  performed  a  pilot-scale  lock  wall  simulation  using  a  fine 
sand  backfill.  The  simulation  was  conceived  to  model  placement  and 
compaction  of  the  backfill  in  a  lock  wall,  surcharge  application,  and 
inundation  of  the  backfill.  SOILSTRUCT-ALPHA  analyses  were  per¬ 
formed  to  model  the  different  stages  of  the  simulation.  Comparisons 
between  SOILSTRUCT-ALPHA  calculations  and  test  measurements 
suggested  that  the  model  can  accurately  predict  the  magnitude  of  the 
downdrag  force  developed  at  the  wall-backfill  interface  in  lock  walls. 

However,  the  applicability  and  accuracy  of  the  model  for  other  types  of 
interfaces  were  not  evaluated.  The  main  objectives  of  this  investigation 
were  the  validation  of  the  model  against  results  of  tests  performed  at  the 
interface  between  coarse  sand  and  concrete,  and  the  expansion  of  the  avail¬ 
able  database  on  interface  hyperbolic  parameters. 

To  accomplish  these  objectives,  it  was  necessary  to  perform  the  follow¬ 
ing  laboratory  activities: 

a.  Select  a  coarse  sand  specimen  for  interface  testing. 

b.  Perform  grain  size  distribution,  minimum/maximum  density,  specific 
gravity,  consolidation,  and  triaxial  testing  on  Blacksburg  Sand. 
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c.  Conduct  initial  loading  tests  to  determine  the  hyperbolic  parameter 
values  of  the  interface  between  concrete  and  dense  Blacksburg  Sand. 

d.  Perform  multidirectional  stress  path  test  to  obtain  data  for  validating 
the  extended  hyperbolic  model. 

A  coarse,  well-graded  sand  (Blacksburg  Sand)  was  selected  for  interface 
testing.  The  basic  properties  of  the  sand  were  determined  from  a  series  of 
laboratory  tests,  such  as  minimum/maximum  density  and  grain  size  analy¬ 
ses.  CD  triaxial  tests  were  also  performed  to  determine  sets  of  hyperbolic 
parameter  values  for  this  soil  (Appendix  A).  These  hyperbolic  parameter 
values  are  consistent  with  values  reported  by  Duncan  et  al.  (1980)  for  soils 
of  similar  gradation. 

Initial  loading  tests  were  performed  at  the  interface  between  concrete 
and  a  coarse,  well-graded  sand.  The  hyperbolic  parameter  values  for  the 
interface  were  determined  based  on  the  results  of  these  tests.  A  multidirec¬ 
tional  stress  path  test  was  performed  on  the  coarse  sand  against  concrete 
interface.  The  purpose  of  this  test  was  to  obtain  data  for  evaluating  the 
accuracy  of  the  extended  hyperbolic  model. 

Comparisons  between  model  calculations  and  the  results  of  the  multi¬ 
directional  stress  path  test  confirmed  that  the  Gomez-Filz-Ebeling  model 
can  predict  accurately  the  response  of  a  variety  of  soil-to-concrete  inter¬ 
faces.  The  principal  advantages  of  the  model  are: 

a.  It  has  a  simple  mathematical  formulation. 

b.  Hyperbolic  parameter  values  for  different  types  of  interfaces  are 
available  in  the  literature. 

c.  It  captures  the  main  features  of  the  interface  response  under  simulta¬ 
neous  changes  in  shear  and  normal  stress  and  unloading-reloading. 

d.  It  provides  accurate  estimates  of  the  response  of  a  variety  of  soil-to- 
concrete  interfaces  subjected  to  complex  stress  paths. 

e.  It  is  relatively  easy  to  implement  in  SSI  analyses. 

/.  It  establishes  a  framework  for  future  work  on  plasticity-based  inter¬ 
face  models. 

The  model  has  several  limitations: 

a.  It  does  not  model  displacements  normal  to  the  interface,  and  the  inter¬ 
face  thickness  is  implicitly  assumed  as  zero. 

b.  It  does  not  model  displacement  softening  of  the  interface. 

c.  The  model  predicts  interface  stiffness  values  that  are  zero  or  negative 
for  certain  loading  combinations.  Although  it  appears  that  negative 
stiffness  values  do  occur  for  certain  loading  combinations,  they  may 
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induce  numerical  problems  in  finite  element  analyses.  Therefore, 
implementation  of  the  model  in  finite  element  programs  requires  the 
use  of  appropriate  stiffness  values  and  numerical  procedures  (see 
Stankowsky,  Runesson,  and  Sture  1993,  for  example)  to  prevent 
these  numerical  problems. 

d.  The  model  has  not  yet  been  used  for  routine  analyses  of  multi- 
anchored  systems.  However,  it  was  used  successfully  by  Gomez, 
Filz,  and  Ebeling  (2000)  to  estimate  vertical  shear  forces  in  a  pilot- 
scale  lock  wall  simulation. 
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Appendix  A 
Results  of  Triaxial  and 
Consolidation  Tests  and 
Determination  of  Hyperbolic 
Parameter  Values 


Triaxial  and  consolidation  tests  were  performed  on  dense  specimens  of 
Blacksburg  Sand.  The  objectives  of  these  tests  were  to 

a.  Determine  the  internal  friction  angle  of  the  sand. 

b.  Determine  the  hyperbolic  parameter  values  of  the  sand. 

c.  Evaluate  hydrocompression  strains  during  inundation. 


A.1  Triaxial  Tests 

The  soil  used  for  all  the  tests  described  herein  is  a  coarse,  well-graded 
sand,  commercially  available  for  the  preparation  of  concrete  in 
Blacksburg,  VA.  This  sand  was  also  used  for  the  interface  tests  described 
in  Chapter  2  and  Appendix  B  of  this  report.  A  more  complete  description 
of  the  properties  of  Blacksburg  Sand  is  presented  in  Chapter  2. 

Drained  triaxial  (CD)  tests  were  performed  to  determine  the  internal 
friction  angle  and  hyperbolic  parameter  values  of  dense  Blacksburg  Sand. 
Three  specimens  were  compacted  to  a  relative  density  of  80  percent.  Each 
specimen  was  subjected  to  an  internal  manometric  pressure  of  -15  to 
-20  kPa,  which  was  gradually  removed  during  application  of  the  cell  pres¬ 
sure.  The  samples  were  de-aired  using  carbon  dioxide,  inundated  with 
de-aired  distilled  water,  and  back-pressure  saturated.  The  specimens  were 
consolidated  under  effective  confining  pressures  ranging  from  69  to 
276  kPa.  Shearing  was  performed  at  a  strain  rate  of  0.25  percent  per 
minute,  which  was  found  to  be  appropriate  for  pore  pressure  dissipation 
during  previous  trials. 
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The  results  of  the  tests  are  presented  graphically  in  Figure  A1 .  All  the 
specimens  exhibited  dilation  during  shear  and  strain-softening  after  mobili¬ 
zation  of  the  peak  strength.  Because  the  peak  strength  envelopes  of  both 
soils  are  curved,  the  value  of  peak  secant  friction  angle,  <j),  for  a  given  con¬ 
fining  stress,  a'3,  can  be  calculated  from  the  following  expression  (Duncan 
etal.  1980): 


4>  =<t>o  —  A<() -logio 


[VI 


(Al) 


where 


<|)o  =  peak  secant  friction  angle  at  a  confining  pressure  of 
101.4  kPa  (1  atm) 

A<(»  =  reduction  in  the  peak  secant  friction  angle  value  for  a  10-fold 
increase  in  a3' 

Determination  of  the  values  of  friction  parameters  <|)o  and  A<|)  is  done 
using  diagrams  of  secant  friction  angle  versus  normalized  confining  pres¬ 
sure  o^'!pa,  as  shown  in  Figure  Ale. 


A.2  Consolidation  Testing 

Consolidation  tests  were  performed  on  dense  specimens  of  Blacksburg 
Sand  to  provide  additional  data  on  their  mechanical  properties  and  deter¬ 
mine  their  susceptibility  to  hydrocompression  (Brandon,  Duncan,  and 
Gardner  1990). 

Two  specimens  of  Blacksburg  Sand  were  prepared  by  compaction  to  a 
relative  density  of  75  percent  in  a  dry  condition.  Each  specimen  was  con¬ 
solidated  under  a  series  of  vertical  stress  increments.  Once  a  predetermined 
stress  was  reached,  the  specimen  was  inundated  after  primary  consolida¬ 
tion  was  attained.  The  results  of  these  tests  are  presented  as  strain  versus 
stress  diagrams  in  Figure  A2. 

During  inundation,  compressive  vertical  strain  was  observed  in  both 
specimens.  The  strain  that  takes  place  during  inundation  includes  deforma¬ 
tion  due  to  secondary  compression  of  the  sand  and  collapse  of  the  structure 
due  to  inundation.  In  retaining  walls  subject  to  changes  in  the  water  eleva¬ 
tion  inside  the  backfill,  hydrocompression-induced  settlements  may  intro¬ 
duce  additional  shear  stresses  at  the  backfill-wall  interface.  Estimates  on 
the  magnitudes  of  hydrocompression-induced  settlements  in  retaining 
walls  can  be  done  as  outlined  in  Appendix  A  of  Gomez,  Filz,  and  Ebeling 
(2000).  The  results  presented  in  Figure  A2  add  to  the  database  on 
hydrocompression  properties  of  soils. 
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a.  Stress-strain  and  volumetric  strain  data 

Figure  A1 .  Results  of  CD  triaxial  tests  on  dense  Blacksburg  Sand  (Sheet  1  of  3) 
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c.  Determination  of  strength  parameters  4>0  and  A<|> 
Figure  A1 .  (Sheet  3  of  3) 
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Figure  A2.  Results  of  consolidation  tests  on  dense  Blacksburg  Sand 
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A.3  Hyperbolic  Parameter  Values  for  Dense 
Blacksburg  Sand 

This  section  contains  a  description  of  the  procedure  suggested  by 
Duncan  et  al.  (1980)  for  the  determination  of  hyperbolic  parameter  values 
for  soils  from  results  of  triaxial  tests.  Hyperbolic  parameter  values  for 
dense  Blacksburg  Sand  are  determined  based  on  the  results  of  the  triaxial 
tests  presented  previously  in  this  appendix. 


A.3.1  Determination  of  hyperbolic  parameter  values 

The  first  step  in  the  determination  of  hyperbolic  parameter  values  is 
checking  for  inconsistencies  in  the  data  from  the  CD  triaxial  tests. 

Figure  A3  shows  the  results  of  the  tests  performed  on  dense  Blacksburg 
Sand.  The  data  points  shown  in  the  figure  are  identical  to  those  shown  in 
Figure  Al.  To  minimize  inconsistencies  in  the  data,  a  smooth  response  of 
the  soil  to  triaxial  testing  was  assumed  for  determination  of  the  hyperbolic 
parameter  values.  This  assumed  response  corresponds  to  the  solid  lines  in 
Figure  A3. 

The  next  step  is  determination  of  the  deviator  stress  at  failure  (a;-C3)y 
for  each  confining  stress.  The  values  of  ( c}-a3)j-  can  be  determined  from 
the  plots  of  deviator  stress  (o;-o3)  versus  axial  strain  e.  Column  2  in  the 
table  presented  in  Figure  A4  contains  the  values  of  deviator  stress  at  fail¬ 
ure  determined  from  Figure  A3.  The  values  of  (aro3)  corresponding  to 
70  and  95  percent  of  (a;-c3)y  are  calculated  as  shown  in  columns  3  and  6, 
respectively. 

The  values  of  axial  strain  e  corresponding  to  70  and  95  percent  of  the 
strength  are  determined  from  the  (0j-03)  versus  e  plots.  Columns  4  and  7 
in  Figure  A4  contain  the  axial  strain  values  determined  as  shown  in 
Figure  A3. 

The  values  in  columns  2,  3,  4,  6,  and  7  are  the  basis  for  the  determina¬ 
tion  of  the  values  of  initial  Young’s  modulus  E(  and  failure  ratio  Rj.  The 
sequence  of  calculations  leading  to  the  determination  of  the  values  of  E( 
and  Rj  is  shown  in  Figure  A4. 

The  value  of  failure  ratio  R ^  to  be  used  for  modeling  is  the  average  of 
the  values  determined  in  Figure  A4.  The  values  of  A  and  n  are  determined 
by  plotting  the  normalized  values  of  initial  Young’s  modulus  against  the 
normalized  confining  stress  in  logarithmic  scale  as  shown  in  Figure  A5. 

Figures  A3  and  A4  also  illustrate  the  determination  of  the  values  of 
bulk  modulus  B  from  the  triaxial  test  data.  If  the  plot  of  volumetric  strain 
ey  versus  axial  strain  e  does  not  reach  a  horizontal  tangent  before  mobiliza¬ 
tion  of  70  percent  of  the  strength,  the  ev  value  corresponding  to  70  percent 
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Volumetric  strain,  sv  (compression  is  positive)  Deviator  stress,  (<Tr< 


Axia  strain,  s 


£=  0.0233  cr'3  =  276  kPa 


0.95  x  (crr(T3)f=  943.4  kPa 


e=  0.0129 


0.7x(  (7?-(73)f  =  695.1  kPa 


o’3  =  138  kPa 


551.6  kPa 


M:  0.0181 

16.4  kPaj 

>6 


//  i 


0.0102 


10.0058 


306.4  kPa 


225.8  kPa 


<73  =  69  kPa 


0  m 
0.0050  T 


0.0040 


£V=  0.00345 


0.0030 


0.00270 


The  volumetric  strain-axial  strain  diagrams 
have  a  horizontal  tangent  before  mobilization 
of  70  percent  of  the  strength.  Therefore,  the 
volumetric  strain  and  deviator  stress  values 
corresponding  to  the  maximum  compression 
of  the  specimens  are  used  for  determination  of 
the  bulk  modulus 


0 

0.0050 


0.0040 


0.0030 


0.0020 


0.0010 


0.00061 


0.0020 


0.0010 


0.0000 


0.0000 


0.00  0.01  0.02  0.03  0.04  0.05 

Axial  strain,  e 


Figure  A3.  Determination  of  axial  and  volumetric  strain  values  at  70  and  95  percent  of  strength.  Data 
from  CD  triaxial  tests  on  dense  Blacksburg  Sand 
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Figure  A4.  Determination  of  the  normalized  values  of  £/  and  B  for  each  of  the  CD  triaxial  tests  performed  on  dense  Blacksburg  Sand  (adapted 
from  Duncan  et  al.  1980) 


Figure  A5.  Determination  of  hyperbolic  parameters  K  and  n  for  dense  Blacksburg  Sand 
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of  the  strength  is  used  for  the  determination  of  B.  If  the  ev  versus  6  plot 
reaches  a  horizontal  tangent  before  mobilization  of  70  percent  of  the 
strength,  the  maximum  value  of  ev  is  used  for  the  determination  of  B. 

The  e  versus  e  plots  of  dense  Light  Castle  Sand  shown  in  Figure  A3 
reach  a  horizontal  tangent  before  mobilization  of  70  percent  of  the 
strength.  The  maximum  values  of  volumetric  strain  ev  are  determined  as 
shown  in  the  figure.  They  are  copied  to  column  10  of  the  table  presented  as 
Figure  A4.  The  deviator  stress  (aro3)  corresponding  to  the  point  of  maxi¬ 
mum  volumetric  strain  is  also  determined  from  the  figure  and  copied  to 
column  9  in  Figure  A4.  These  two  values  are  used  to  determine  the  value 
of  B  for  each  of  the  confining  stresses  applied.  The  values  of  Kb  and  m  are 
determined  from  a  logarithmic  plot  of  normalized  bulk  modulus  versus 
normalized  confining  stress,  as  shown  in  Figure  A 6. 

Table  A1  summarizes  the  values  of  the  hyperbolic  parameters  K,  n,  Rp 
Kb,  and  m  for  dense  Blacksburg  Sand,  as  well  as  the  values  of  and  A<t> 
determined  from  the  triaxial  tests  data,  as  described  previously  in  this 
appendix.  Parameter  values  obtained  from  the  results  of  drained  CD 
triaxial  tests  on  a  variety  of  soils  are  reported  in  Table  5  of  Duncan  et  al. 
(1980).  These  values  are  reproduced  in  Section  A. 3. 3. 


Table  Al 

Hyperbolic  Parameter  Values  for  Dense  Blacksburg  Sand 
(Dr  =80%) 

Parameter1 

Value 

K 

642 

n 

0.25 

R, _ _ _ 

0.55 

Kh  _ 

567 

— - — — — ■ — 

m 

0.01 

<t>„.  deg 

43.4 

Adi,  deg 

7.4 

1  Parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 

As  discussed  in  Chapter  2  of  this  report,  the  values  of  the  hyperbolic 
parameters  for  dense  Blacksburg  Sand  appear  consistent  with  those 
reported  by  Duncan  et  al.  (1980). 


A.3.2  Comparison  of  model  to  test  data 

The  stress-strain  response  from  the  model  is  calculated  using  the 
following  expression: 
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Figure  A6.  Determination  of  hyperbolic  parameters  Kb  and  m  for  dense  Blacksburg  Sand 
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(oi-®3)=- 


K  Pa' 


(„'\n 

<Paj 


-  +  R 


/- 


(O1-O3) 


(A2) 


The  deviator  stress  at  failure  ( ora3 )y  is  calculated  from  the  following 
expression: 


/  x  o3  sm$ 

(0,  -03 ),  =2— - - 

'  'f  1  -  sm(J) 


(A3) 


where  <j)  can  be  determined  from  Equation  A1 . 

The  volumetric  strain-axial  strain  response  is  calculated  using  the 
following  expression: 


ey  = 


3 -Kbpa- 


(„'\m 

q3 

<Pa. 


(A4) 


The  hyperbolic  stress-strain  and  volumetric  strain-axial  strain  response 
of  dense  Blacksburg  Sand  was  calculated  using  Equations  A2,  A3,  and  A4 
and  the  hyperbolic  parameters  determined  following  the  procedure 
described  previously.  Figure  A7  shows  a  comparison  between  the  test  data 
and  the  calculated  hyperbolic  response.  In  the  figure,  the  stress-strain 
hyperbolas  are  interrupted  at  the  value  of  deviator  stress  at  failure 
A  horizontal  stress-strain  relationship  (i.e.,  zero  Young’s  modulus)  is  used 
to  model  the  response  of  the  soil  at  failure. 

It  can  be  seen  that  the  hyperbolic  model  provides  an  accurate  approxi¬ 
mation  of  the  stress-strain  response  measured  during  the  tests.  The  volu¬ 
metric  strain-axial  strain  response  calculated  using  the  hyperbolic  model 
also  provides  a  good  approximation  of  the  test  data  for  the  initial  stages  of 
shear,  in  which  compression  takes  place.  It  does  not  model  subsequent 
dilation  of  the  soil. 
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Volumetric  strain,  ev  (compression  is  positive)  Deviator  stress,  (aro3)  (kPa) 


Axial  strain,  s 


0.00  0.01  0.02  0.03  0.04  0.05 


0.00  0.01  0.02  0.03  0.04  0.05 

Axial  strain,  e 


Figure  A7.  Hyperbolic  model  for  dense  Blacksburg  Sand  and  comparison  to  CD  test  data 
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A.3.3  Selection  of  hyperbolic  parameter  values  for  soils  under 
drained  conditions 

A  summary  of  property  values  for  different  types  of  soils  tested  under 
drained  conditions  is  presented  by  Duncan  et  al.  (1980).  The  summary  is 
reproduced  in  Table  A2.  Additionally,  Stark,  Ebeling,  and  Vettel  (1994) 
provide  hyperbolic  stress-strain  parameters  for  freshly  deposited,  normally 
consolidated  silts  and  clayey  silts,  which  are  reproduced  in  Tables  A3  and 
A4.  Note  that  the  unit  weight  given  in  Tables  A3  and  A4  for  the 
kaolinite-silt  and  montmorillonite-silt  specimens  is  the  dry  unit  weight, 

Ydry.  The  moist  unit  weight,  ymoist,  is  computed  using  ymoist  =  ydry  (1  +  w) 
where  the  water  content,  w,  is  specified  as  a  fraction. 

The  data  contained  in  this  section  are  not  intended  to  replace  site- 
specific  testing  on  high-quality  soil  specimens,  but  are  for  use  as  interim 
guidance  prior  to  the  availability  of  laboratory  test  results.  The  data  in 
Tables  A2  through  A4  can  be  used  as  an  aid  for  selecting  hyperbolic 
parameter  values  for  analysis  when  test  results  are  not  available. 

Average  hyperbolic  parameter  values  for  coarse  soils  were  determined 
from  the  data  in  Table  A2  and  are  presented  in  Table  A5.  The  table  also 
indicates  the  number  of  data  points  from  Table  A2  that  were  used  for  the 
determination  of  the  average  values  for  each  soil  type.  To  provide  an  idea 
of  the  variability  of  hyperbolic  parameter  values,  the  standard  deviation  is 
also  indicated  in  parentheses  after  the  average  value  for  each  of  the 
parameters. 

Based  on  the  data  reproduced  in  Table  A2,  Duncan  et  al.  (1980)  pre¬ 
sented  a  conservative  interpretation  of  parameter  values  for  various  types 
of  soils  at  different  degrees  of  compaction  relative  to  Standard  Proctor. 
These  values  are  reproduced  in  Table  A6.  The  original  1980  version  of  the 
table  provides  modulus  and  strength  values  that  are  relatively  low,  and  unit 
weights  (single  value)  that  are  relatively  high  for  the  types  of  materials 
listed.  Consequently,  the  use  of  these  values  in  finite  element  analyses  may 
tend  to  provide  conservative  (large)  deformation  values.  This  table  was 
updated  by  Idriss  and  Duncan  (1988)  to  include  ranges  of  unit  weights,  and 
these  ranges  are  also  listed  in  Table  A6. 

Table  A7  provides  ranges  of  hyperbolic  parameter  values  for  coarse 
soils.  These  ranges  are  based  on  the  data  presented  in  Table  A2,  combined 
with  the  judgment  and  experience  of  the  authors  of  this  report.  In  the 
interim  and  until  data  from  laboratory  tests  are  available,  Table  A7  may  be 
used  as  a  starting  point  for  selecting  soil  property  values  for  soil-structure 
interaction  analyses.  However,  it  is  recommended  that  the  analyst  also 
refer  directly  to  Table  A2  when  selecting  property  values  because  factors 
such  as  mineralogy,  angularity,  etc.,  which  are  not  reflected  in  Table  A7, 
can  influence  soil  behavior.  Further,  the  purposes  of  the  analyses  and  the 
nature  of  the  applied  loading  should  be  considered  when  selecting  property 
values.  One  of  the  best  procedures,  when  feasible,  is  to  calibrate  property 
values  and  analysis  procedures  with  observed  performance. 
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Table  A2 

Stress-Strain  and  Strength  Parameters  for  Soils  Tested  Under  Drained  Conditions  (After 
Table  5,  Duncan  et  al.  1980) _ 


Siltv  Sandy  Gravel  (Oroville  Dam)  T9 


Amphibolite  Gravel  (Oroville  Dam  Shell)  [14 


Crushed  Basaltic  Rock  (Round  Butte  Dam)  [17 


Sandy  Gravel  (Rowallan  Dam)  [4 


Clayey  Gravel  (New  Hogan  Dam  Core)  [2 


Argillite  Rockfill  (Pyramid  Dam  Shell)  [14 


Crushed  Olivine  Basalt  [14 


Siltv  Sand,  Some  Gravel  (Round  Butte  Dam)  [17 


Venato  Sandstone  (0.5  in.  max.  size)  [1 


Glacial  Cutwash  Sand  [10 


Sacramento  River  Sand  [13 


Sacramento  River  Sand  [13 


Sacramento  River  Sand  [1 3 


Sacramento  River  Sand  [13 


Poorly  Graded  Sand  (Port  Allen  Lock)  [20 


Poorly  Graded  Sand  (Port  Allen  Lock)  [20 


Poorly  Graded  Sand  (Port  Allen  Lock)  [20 


4.1 

1.8 

0.6 

1.7 

0.09 

0.009 

0.17 

0.07 

0.025 

0.03 

0.4 

0.14 

0.22 

0.17 

0.15 

0.22 

0.17 

0.15 

0.22 

0.17 

0.15 

0,22 

0.17 

0.15 

0.25 

0.17 

0.15 

0.25 

0.17 

0.15 

0.2 

0.17 

0.12 

0.2 

0.17 

0.12 

0.2 

0.17 

0.12 

(Sheet  1  of  7) 


Reference  citations  are  provided  at  the  conclusion  of  Table  A2. 
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Table  A2  (Continued 


Compaction 


Row 

No.  Grou 


GW-1 


GW-2 


GW-3 


GW-4 


GW-5 


GW-6 


GW-7 


GP-2 


GP-3 


GP-6 


GP-7 


GP-11 


GP-13 


GC-1 


SW-1 


SW-2 


SW-3 


SW-5 


SP-3 


SP-4A 


SP-4B 


SP-4C 


SP-4D 


SP-5A 


SP-5B 


SP-7A 


SP-7B 


SP-7C 


16,450  120.0 


118.3 


100.0 

13  0 

95.5 

100.0 

13.0 

100.0 

100.0 

13.0 

105.1 

Init. 

Void 

Ratio 

Relative 

Density(%) 

0.39 

70 

0.32 

95 

133.8 


148.0 


152.0 


99.0 


135.0 


10.8  107.0  10.8 


111.6 


125.4 


13.2  108.7  13.5 


117.5 


112.3 


89.5 


94.0 


97.8 


103.9 


65 


50 


50 


95 


100 


100 


99 


0  233  100 


0.46 

100 

0.43 

100 

0.47 

93 

0.5 

80 

0.87 

38 

IEH 

60 

0.71 

78 

0.61 

100 

0.82 

Loose 

0.64 

Dense 

0.73 

49 

0.65 

73 

0,57 

98 

Sub-anqular 


Sub-rounded 


Sub-rounded 


Sub-rounded 


Sub-rounded 


Sub-anqular 


Anqular 


Rounded 


Rounded 


Anqular 


Rounded 


Anqular 


Sub-rounded 


Sub-rounded 


Rounded 


Rounded 


Rounded 


Rounded 


Rounded 


Rounded 


Rounded 


Rounded 


Rounded 


(Sheet  2  of  7) 


'Ratings  were  based  on  estimated  consistency  and  reliability  of  data,  judged  by  comparing  stress-strain  and  volume  change  curves 
calculated  using  the  listed  parameters  with  the  experimental  curves.  Definitions  of  ratings  are  as  follows:  ***  =  excellent,  -  very 
good,  *  =  good.  Data  of  lesser  quality  were  discarded. 
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Table  A2  (Continued) 

Row 

No. 

Soil 

Soil  Description  [Referencel1 

Grain  Size  (mm) 

LL 

PI 

HS 

D10 

29 

SP 

SP-12 

Coarse  to  Fine  Sand  (Round  Butte  Dam)  [171 

minus  No.  4  s 

eve 

NP 

NP 

30 

SP 

SP-13 

Pumicecus  Sand  (Round  Butte  Dam)  [171 

0.85 

0.41 

0.24 

| 

31 

SP 

SP-14 

Pumicecus  Sand  (Round  Butte  Dam)  [191 

1.0 

0.5 

0.24 

| 

1 

32 

SP 

SP-16A 

Fine  Silica  Sand  (Loose)  [81 

0.27 

0.2 

0.165 

33 

SP 

SP-16B 

Fine  Silica  Sand  (Dense)  [81 

0.27 

0.2 

0.165 

34 

SP 

SP-17A 

Monterey  No.  0  Sand  (Cylind.  specimen)  [121 

0.43 

0.37 

0.29 

NP 

NP 

35 

SP 

mm 

Monterev  No.  0  Sand  (Cubical  specimen)  [121 

0.43 

0.37 

0.29 

NP 

NP 

36 

SP 

SP-17C 

Monterev  No.  0  Sand  (Cylind.  specimen)  [121 

0.43 

0.37 

0.29 

NP 

NP 

37 

SP 

SP-17D 

Monterev  No.  0  Sand  (Cubical  specimen)  [121 

0.43 

0.37 

0.29 

NP 

NP 

38 

SP 

SP-18 

Basaltic  Sand  (Round  Butte  Dam)  [181 

3. 

9. 

0.13 

39 

SM 

SM-4 

Siltv  Sand  (Chatfield  Dam)  [211 

0.62 

0.16 

0.026 

20 

0 

40 

SM 

SM-5 

Silty  Gravelly  Sand  (Chatfield  Dam)  [211 

1.15 

0.28 

0.05 

NP 

NP 

41 

SM 

SM-6 

Siltv  Sand  w/Pebbles  (Round  Butte  Dam)  1171 

0.31 

0.1 

0.04 

NP 

NP 

42 

SM 

SM-9 

Siltv  Sand  w/Pumice  (Round  Butte  Dam)  [17] 

0.15 

0.054 

0.013 

NP 

NP 

43 

SM 

SM-13 

Siltv  Sand  (Round  Butte  Dam)  [191 

0.27 

0.027 

0.0022 

44 

SM 

SM-16 

Siltv  Sand  and  Gravel  (Round  Butte  Dam)  [181 

0.45 

0.052 

45 

SM-SC 

SM-SC-1A 

Siltv  Clayey  Sand  (Mica  Dam  Core)  [5,  111 

0.34 

0.03 

0.002 

21 

HI 

46 

SM-SC 

SM-SC-1 B 

Siltv  Clayey  Sand  (Mica  Dam  Core) [5,  111 

0.34 

0,03 

0.002 

21 

11 

47 

SM-SC 

SM-SC-1 C 

Silty  Clayey  Sand  (Mica  Dam  Core)  [5,  111 

0.34 

0.03 

0.0002 

21 

H 

48 

ML 

ML-1 

Cannonsville  Silt  (Undisturbed)  [101 

0.033 

0.018 

0.005 

49 

ML 

ML-4 

Sandy  Silty  w/Pumice  (Round  Butte  Dam)  [171 

0.078 

0.032 

0.0064 

NP 

NP 

50 

ML 

ML-5 

Sandy  Silty  w/Pumice  (Round  Butte  Dam)  [171 

0.1 

0.025 

0,0052 

NP 

NP 

51 

CL 

CL-29C 

Siltv  Clay  (Canyon  Dam)  [71 

0.037 

0.008 

— 

34 

19 

52 

CL 

CL-29D 

Siltv  Clay  (Canyon  Dam)  [71 

0.037 

0.008 

— 

34 

19 

53 

CL 

CL-30E 

Silty  Clay  (Canyon  Dam)  [71 

0.037 

0.008 

— 

34 

19 

54 

CL 

CL-30F 

Siltv  Clay  (Canyon  Dam)  [7] 

0.037 

0.008 

- 

34 

19 

55 

CL 

CL-34E 

Siltv  Clay  (Canyon  Darn)  [61 

0.037 

0.008 

- 

34 

19 
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1  Reference  citations  are  provided  at  the  conclusion  of  Table  A2. 


Appendix  A  Results  of  Triaxial  and  Consolidation  Tests  and  Determination  of  Hyperbolic  Parameter  Values 


A19 


Table  A2  (Continued 


29 

SP-12 

30 

SP-13 

31 

SP-14 

32 

SP-16A 

33 

SP-16B 

34 

SP-17A 

35 

SP-17B 

36 

SP-17C 

37 

SP-17D 

38 

SP-18 

39 

SM-4 

40 

SM-5 

41 

SM-6 

42 

SM-9 

43 

SM-13 

44 

SM-16 

45 

SM-SC-1  A 

46 

SM-SC-1B 

47 

SM-SC-1  C 

48 

ML-1 

49 

ML-4 

Init.  Relative 

w/c  Void  Density 

(%)  Ratio  (%) 


120.1 

9.5 

120.0 

123.0 

9.5 

116.7 

132  0 

8.1 

124.5 

110  6 

17.5 

108.1 

91 .7 

19.5 

88.4 

105  6 

16.4 

104.5 

109.3 

12.9 

109.0 

Std.AASHO  136.0 


Std.AASHO  136.0 


Std.AASHO  136.0 


16,450 


16,450 


Harvard 


Harvard 


Harvard 


Harvard 


Harvard 


134.0 


128.2  11.9 


108.0 


92.8  17.7 


99.2  17.0 


111.2  13,1 


116.2  13.3 


115.1  15.2 


110.0  17.4 


106.3  19.0 


(Sheet  5  of  7) 


’Ratings  were  based  on  estimated  consistency  and  reliability  of  data,  judged  by  comparing  stress-strain  and  volume  change  curves 
calculated  using  the  listed  parameters  with  the  experimental  curves.  Definitions  of  ratings  are  as  follows:  ***  =  excellent,  **  =  very 
good,  *  =  good.  Data  of  lesser  quality  were  discarded. 
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Table  A3 

Drained  Hyperbolic  Stress-Strain  Parameters  for  Kaolinite-Silt  Mixtures 
Consolidated  Drained  Triaxial  Tests  Listed  in  Table  3  of  Stark,  Ebeling, 

from  Is 
and  Vet 

otropic: 
tel  (199 

ally 

4) 

Percent 

clay 

mineral9 

(1) 

Standard 

Proctor 

relative 

compaction 

(%) 

(2) 

Average 
initial 
dry  unit 
weight 
(kN/m3) 

(3) 

Average 

Initial 

water 

content 

(percent) 

(4) 

Range  of 
triaxial 

consolidation 

pressure 

(kPa) 

(5) 

Effective 

stress 

cohesion 

(KPa) 

(6) 

Effective 

stress 

friction 

angle 

(deg) 

(7) 

Modulus 

number 

K 

(8) 

Modulus 

exponent 

n 

(9) 

Bulk 

modulus 

number 

xb 

(10) 

Bulk 

modulus 

exponent 

m 

(11) 

Failure 

ratio 

(12) 

mm 

m 

1,120-1,740 

0 

■ 

l  1 

mm 

755-1,725 

0 

HI 

1.00 

0%  Kao 

28 

765-1 ,585 

0 

35 

9 

H 

1.00 

0%  Kao 

85 

N/A 

N/A 

N/A 

1^1 

H 

10%  Kao 

0 

37 

10%  Kao 

0 

10%  Kao 

90 

23 

0 

10%  Kao 

85 

300-1,390 

0 

30%  Kao 

100 

ESI 

17 

980-1,685 

0 

33 

105 

1.00 

9 

30%  Kao 

95 

■ 

795-1,485 

31 

70 

H 

H 

1.00 

mill 

30%  Kao 

90 

■ 

395-1,520 

0 

30 

m 

■ 

H 

mom 

30%  Kao 

85 

17.9 

m 

300-1 ,660 

0 

29 

60 

H 

H 

50%  Kao 

17.4 

18 

880-1,670 

0 

28 

65 

1.0 

30 

1.00 

9 

50%  Kao 

17.1 

19 

785-1,475 

0 

27 

60 

1.0 

25 

1.00 

0.70 

50%  Kao 

mm 

20 

390-1,275 

0 

26 

55 

1.0 

20 

1.00 

0.65 

50%  Kao 

mm 

21 

350-1,280 

0 

25 

50 

1.0 

15 

1.00 

0.60 

aKao  =  kc 
bN/A  =  N< 

lolinite. 
at  available. 

Appendix  A  Results  of  Triaxial  and  Consolidation  Tests  and  Determination  of  Hyperbolic  Parameter  Values 


A23 


Table  A4 

Drained  Hyperbolic  Stress-Strain  Parameters  for  Montmorillonite-Silt  Mixtures  from 
Isotropically  Consolidated  Drained  Triaxial  Tests  Listed  in  Table  5  of  Stark,  Ebeling,  and 
Vettel  (1994) 


Standard 

Average 

Average 

Range  of 

Effective 

Proctor 

initial 

initial 

triaxial 

Effective 

stress 

Bulk 

Bulk 

Percent 

relative 

dry  unit 

water 

consolidation 

stress 

friction 

Modulus 

Modulus 

modulus 

modulus 

Failure 

clay 

compaction 

weight 

content 

pressure 

cohesion 

angle 

number 

exponent 

number 

exponent 

ratio 

mineral3 

(%) 

(kN/m3) 

(%) 

(kPa) 

(kPa) 

(deg) 

K 

n 

Kb 

m 

*, 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

15.4  27  200-1,640  0  40  270  1.0 


Mont  =  montmorillonite. 


16.0  24  1,080-1,620  0 


15.7  26  990-1,620  0  20 


90  1.0 


1.0  0.75 


Table  A5 

Average  Property  Values  of  Coarse  and  Fine  Soils  in  Drained  Condition 


Notes: 

Symbols  are  defined  in  the  notation  (Appendix  C). 
Numbers  in  parentheses  indicate  the  standard  deviation. 
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Table  A6 

Property  Values  for  Compacted  Soils  Listed  in  Table  7  of  Duncan  et  al.  (1980) 


Range1 


Unified  Dr 

Classification  (Vo)1 


RC 

<%) 


(kN/m3) 


in?m 

(kN/m3) 


<t> 

(deg) 


A<(> 

(deg) 


(kPa) 


K 


n 


K. 


m 


GW,  GP,  SW,  SP 


100 

75 

50 

25 


105 

100 

95 

90 


23.6 

22.8 

22.0 

21.2 


16.5-23.6 

15.7-22.8 

14.9-22.0 

14.1-21.2 


42 

39 

36 

33 


9 

7 

5 

3 


0 

0 

0 

0 


600 

450 

300 

200 


0.4 

0.4 

0.4 

0.4 


0.7 

0.7 

0.7 

0.7 


175 

125 

75 

50 


0.2 

0.2 

0.2 

0.2 


SM 


100 

95 

90 

85 


21.2 

20.4 

19.6 

18.9 


17.3-21.2 

16.5-20.4 

15.7-19.6 

14.9-18.9 


36 

34 

32 

30 


8 

6 

4 

2 


0 

0 

0 

0 


600 

450 

300 

150 


0.25 

0.25 

0.25 

0.25 


0.7 

0.7 

0.7 

0.7 


450 

350 

250 

150 


0.0 

0.0 

0.0 

0.0 


SM-SC 


100 

95 

90 

85 


21.2 

20.4 

19.6 

18.9 


18.1-21.2 

17.3-20.4 

16.5-19.6 

15.7-18.9 


33 

33 

33 

33 


0 

0 

0 

0 


0.5 

400 

0.6 

0.4 

200 

0.6 

0.3 

150 

0.6 

0.2 

100 

0.6 

0.7 

200 

0.5 

0.7 

100 

0.5 

0.7 

75 

0.5 

0.7 

50 

0.5 

CL 


100 

95 

90 

85 


21.2 

20.4 

19.6 

18.9 


18.1-21.2 

17.3-20.4 

16.5-19.6 

15.7-18.9 


30 

30 

30 

30 


0 

0 

0 

0 


0.4 

0.3 

0.2 

0.1 


150 

120 

90 

60 


0.45 

0.45 

0.45 

0.45 


0.7 

0.7 

0.7 

0.7 


140 

110 

80 

50 


0.2 

0.2 

0.2 

0.2 


Note: 


Parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 

RC  refers  to  the  relative  compaction  according  to  Standard  Proctor:  ASTM  D698-78. 


ym  refers  to  moist  unit  weight. 

1  Listed  in  Table  8-1,  Chapter  8  (Idriss  and  Duncan  1988). 


Table  A7 

Proposed  Ranges  of  Starting  Property  Values  for  Granular  Soils  in  Drained  Condition 


Unified 

Classification 

ft) 

(deg) 

A<|> 

(deg) 

(kPa) 

K 

n 

m 

04 

0 

200-600 

0.4-0.6 

0.7-0.9 

150-300 

0. 1-0.4 

0 

l 

0.7-0. 9 

200-400 

0. 1-0.4 

40-50 

mm 

0 

0.7-0. 9 

400-1100 

0. 1-0.4 

Note:  Symbols  are  defined  in  the  notation  (Appendix  C). 
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In  SOILSTRUCT-ALPHA,  the  soil  property  values  for  primary  loading 
required  as  input  for  analysis  are  the  modulus  number  K,  the  modulus 
exponent  n,  the  failure  ratio  R \p  the  friction  angle  <)),  and  the  nominal  Pois¬ 
son’s  ratio  vnom.  Starting  values  of  K,  n,  and  R,  can  be  determined  using 
Table  A7,  and  may  be  adjusted  depending  on  the  characteristics  of  the 
problem  under  analysis  and  the  expected  response  of  the  soil.  The  value  of 
the  friction  angle  can  be  determined  using  Equation  A1  and  the  values  of 
<|)0  and  A<j>  from  Table  A7,  together  with  a  representative  value  of  confining 
stress. 


The  value  of  v  can  be  related  to  the  bulk  modulus  number  K,  and 

norm  o 

modulus  number  K,  following  the  procedure  described  in  Appendix  C  of 
Ebeling,  Pace,  and  Morrison  (1997).  This  procedure  provides  a  good  start¬ 
ing  value  of  .V  that  can  be  refined  using  the  one-dimensional  (1-D) 
column  analysis  procedure  as  described  by  Ebeling  and  Wahl  (1997). 

Before  performing  a  1-D  column  analysis,  it  is  necessary  to  estimate  the 
value  of  the  coefficient  of  earth  pressure  at  rest  KQ.  For  a  tall  backfill  of 
homogeneous  granular  fill  where  the  overburden  pressures  exceed  the 
compaction-induced  preconsolidation  pressures  over  the  majority  of  the 
backfill,  an  adequate  value  of  KQ  can  be  estimated  using  Jaky’s  relation¬ 
ship  ( Kg  =  1  -  sin  <t>').  In  relatively  short  backfills  or  in  backfills  compacted 
with  heavy  equipment,  the  overburden  pressures  may  not  exceed  the 
compaction-induced  preconsolidation  pressures  over  the  majority  of  the 
backfill  height  and  the  value  of  K0  may  vary  significantly  over  the  height 
of  the  backfill.  In  these  cases,  the  value  of  Kg  for  analysis  must  account  for 
overconsolidation  of  the  backfill. 


The  value  of  vwom  is  dependent  on  Ka  and  the  stress  level  SL.  In  a  1-D 
column  analysis,  the  value  of  vnom  is  adjusted  until  matching  the  target  KQ 
value  determined  from  Jaky’s  relationship.  A  good  initial  approximation  to 
the  value  of  vnom  can  be  obtained  solving  Equation  A1  from  Appendix  A 
of  Ebeling  and  Wahl  (1997),  which  is  restated  as  Equation  A5  below. 


l-0-2v„om) 


K„  =- 


1-R 


(l-sinM-*,) 
2 Ka  sin<j) 


-|2 


1+0-2V-) 


1  -R 


(l-sin((>)(l-^0) 
2 Ka  sin<|> 


(A5) 


For  the  McAlpine  lock  wall  analyses  described  by  Ebeling  and  Wahl 
(1997),  a  value  KQ  of  0.44  was  estimated  using  Jaky’s  relationship  with 
<t>'  =  34  deg.  The  value  of  failure  ratio  R^oi  the  backfill  was  0.7.  Using 
these  values  and  solving  Equation  A1  in  Ebeling  and  Wahl  (1997)  yielded 
a  value  of  v nom  of  0.025.  In  addition,  1-D  column  analyses  of  the  backfill 
were  performed  using  the  backfill  placement  option  incorporated  in 
SOILSTRUCT-ALPHA.  The  value  of  vnom  obtained  from  the  1  -D  column 
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analyses  was  consistent  with  the  value  of  0.025  determined  from 
Equation  A5. 


The  value  of  Poisson’s  ratio  v  corresponding  to  a  given  value  of  nomi¬ 
nal  Poisson’s  ratio  v/wm  can  be  determined  using  Equation  A2  in  Ebeling 
and  Wahl  (1997).  According  to  this  equation  (restated  as  Equation  A6 
below),  the  value  of  v  depends  on  vnom,  Rf  and  the  stress  level  SL.  For 
v  =  0  025  Rr  =  0.7,  and  SL  -  0.51,  the  value  of  Poisson’s  ratio  v  is  0.3. 

nom  ’  J  ’ 


(A6) 


For  unloading-reloading,  SOILSTRUCT-ALPHA  requires  an  additional 
property  value,  the  unload-reload  modulus  Kur.  Values  of  the 
unload-reload  modulus  number  Kur  for  soils  are  not  readily  available  in  the 
literature.  According  to  Duncan  et  al.  (1980),  the  value  of  Kur  is  always 
larger  than  the  value  of  the  modulus  number  K  for  primary  loading.  For 
stiff  soils  such  as  dense  sands,  Kur  may  be  20  percent  greater  than  K.  For 
soft  soils  such  as  loose  sands,  Kur  may  be  three  times  as  large  as  K. 

A  more  precise  determination  of  the  value  of  Kur  can  be  achieved  from 
the  results  of  a  1-D  consolidation  test  that  includes  one  or  more  cycles  of 
unloading-reloading.  A  procedure  for  the  determination  of  Kur  from  1-D 
unloading-reloading  tests  is  presented  in  Chapter  5  of  Clough  and  Duncan 
(1969). 
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Appendix  B 

Results  of  Interface  Tests  and 
Determination  of  interface 
Hyperbolic  Parameter  Values 


This  appendix  contains  the  results  of  the  interface  tests  performed  for 
this  investigation.  It  also  contains  a  brief  description  of  the  procedure  de¬ 
veloped  by  Clough  and  Duncan  (1971)  for  the  determination  of  hyperbolic 
parameter  values  for  interfaces.  Hyperbolic  parameter  values  of  the  inter¬ 
face  between  concrete  and  dense  Blacksburg  Sand  are  determined  based  on 
the  Clough  and  Duncan  (1971)  procedure. 


B.1  Results  of  Interface  Tests 

The  results  of  the  interface  tests  performed  for  this  investigation  are  pre¬ 
sented  in  Figures  B1  through  B3.  The  results  of  initial  loading  tests 
T2001_5,  T2002_10,  T2003_20,  and  T2004_40  are  presented  in  Figures  B1 
and  B2.  During  these  tests,  the  interface  was  sheared  under  constant  nor¬ 
mal  stress  until  failure  was  attained.  Determination  of  the  interface  hyper¬ 
bolic  parameter  values,  which  is  described  in  the  following  section,  was 
done  based  on  the  results  of  these  tests. 

Figure  B3  shows  the  results  of  multidirectional  stress  path  test  T2005. 
The  test  consisted  of  the  application  of  a  stress  path  that  included  simulta¬ 
neous  changes  in  shear  and  normal  stresses  and  unloading-reloading.  The 
purpose  of  this  test  was  to  provide  a  basis  for  a  performance  evaluation  of 
the  extended  hyperbolic  model  under  complicated  loading  paths. 
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Interface  shear  stress,  x  (kPa) 
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Interface  shear  stress,  t  (kPa)  I  ■  I  Interface  shear  stress,  x  (kPa) 


0  25  50  75  100  125  1 
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b.  Shear  stress  versus  interface  displacement  data 

Figure  B3.  Multidirectional  stress  path  test  T2005  on  dense  Blacksburg 
Sand-to-concrete  interface 
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B.2  Determination  of  Interface  Hyperbolic 
Parameter  Values 

The  first  step  in  the  determination  of  hyperbolic  parameter  values  is 
checking  for  inconsistencies  in  the  data  from  the  interface  tests.  Figure  B4 
shows  the  results  of  the  tests  performed  on  the  dense  Blacksburg  Sand 
against  concrete  interface.  The  data  points  shown  in  the  figure  are  identical 
to  those  shown  in  Figure  B 1 .  Frequently,  data  from  interface  tests  may  pre¬ 
sent  some  inconsistencies  that  can  be  minimized  by  developing  an  assumed 
shear  stress-displacement  response.  The  response  assumed  for  the  determi¬ 
nation  of  the  hyperbolic  parameters  in  this  example  is  represented  as  solid 
lines  in  Figure  B4.  Because  of  the  good  precision  of  displacement  and 
shear  stress  measurements  in  these  tests,  few  inconsistencies  can  be  seen  in 
the  data,  and  the  assumed  response  is  practically  identical  to  the  measured 
response. 

The  next  step  is  the  determination  of  the  shear  stress  at  failure  Xyfor 
each  normal  stress.  The  values  of  Xycan  be  determined  from  the  shear 
stress-displacement  plots  of  the  tests.  Column  2  of  the  table  presented  as 
Figure  B5  contains  the  values  of  shear  stress  at  failure  determined  from 
Figure  B4.  The  values  of  shear  stress  corresponding  to  70  and  95  percent 
of  Xy  are  calculated  as  shown  in  columns  3  and  6,  respectively. 

The  values  of  interface  displacement  corresponding  to  70  and  95  per¬ 
cent  of  the  shear  stress  at  failure  are  determined  from  the  shear  stress- 
displacement  plots.  Columns  4  and  7  in  Figure  B5  contain  the 
displacement  values  determined  as  shown  in  Figure  B4. 

The  values  in  columns  2,  3,  4,  6,  and  7  are  the  basis  for  the  determina¬ 
tion  of  the  values  of  initial  interface  stiffness  Ks-  and  failure  ratio  R^.  The 
sequence  of  calculations  leading  to  the  determination  of  the  values  of  Ksi 
and  Rj-j  is  shown  in  Figure  B5. 

The  value  of  failure  ratio  R^  to  be  used  for  modeling  is  the  average  of 
the  values  determined  in  Figure  B5.  The  values  of  Kj  and  tij  are  determined 
by  plotting  the  normalized  values  of  initial  interface  stiffness  against  the 
normalized  normal  stress  in  logarithmic  scale,  as  shown  in  Figure  B6. 

Values  of  interface  hyperbolic  parameters,  determined  following  the  pro¬ 
cedure  outlined  in  this  section,  may  require  further  adjustments  when  using 
the  Gomez-Filz-Ebeling  interface  model  (G6mez,  Filz,  and  Ebeling  2000). 
As  discussed  in  a  subsequent  section,  only  the  value  of  Kj  for  the  interface 
between  dense  Blacksburg  Sand  and  concrete  required  further  adjustment. 
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Interface  shear  stress,  x  (kPa) 


Figure  B4.  Example  determination  of  interface  displacements  at  70  and 
95  percent  of  strength  data  from  interface  tests  on  dense 
Blacksburg  Sand 
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B.3  Comparison  of  Model  to  Test  Data 


The  shear  stress-displacement  response  from  the  model  is  calculated  us¬ 
ing  the  following  expression: 


X  = 


A , 


_ 1 _  | 

(a  Y;  a„  tan8 

Kryw-  — 

\Pa) 


(Bl) 


where  8  is  the  interface  friction  angle.  This  expression  is  valid  for  fric¬ 
tional  interfaces  with  zero  adhesion  intercept. 

The  shear  stress-displacement  response  of  the  interface  between  dense 
Blacksburg  Sand  and  concrete  was  calculated  using  Equation  Bl  and  the 
hyperbolic  parameters  determined  as  described  previously.  Figure  B7  com¬ 
pares  the  test  data  and  the  calculated  hyperbolic  response.  In  the  figure,  the 
shear  stress-displacement  hyperbolas  are  interrupted  at  the  value  of  shear 
stress  at  failure,  Xr.  A  horizontal  shear  stress-displacement  relationship 
(i.e.,  zero  interface  stiffness)  is  used  to  model  the  response  of  the  soil  after 
failure  is  attained.  It  can  be  seen  that  the  hyperbolic  model  provides  an  ac¬ 
curate  approximation  of  the  interface  response  measured  during  each  of  the 
tests. 


B.4  Adjustment  of  the  Hyperbolic  Parameter 
Values 

As  observed  by  Gomez,  Filz,  and  Ebeling  (2000),  hyperbolic  parameter 
values,  determined  according  to  the  method  outlined  in  the  previous  sec¬ 
tions,  are  usually  accurate  enough  for  implementation  of  the  Clough  and 
Duncan  (1971)  hyperbolic  model  for  interface  shear  under  constant  normal 
stress.  However,  for  implementation  of  the  G6mez-Filz-Ebeling  model  for 
interfaces,  it  may  be  necessary  to  improve  the  accuracy  of  the  parameter 
values  through  an  adjustment  procedure. 

The  procedure  for  adjustment  of  hyperbolic  parameter  values  is  de¬ 
scribed  in  detail  by  G6mez,  Filz,  and  Ebeling  (2000).  It  is  based  on  the  nor¬ 
malization  of  the  measured  and  calculated  response  of  the  interface  during 
initial  loading.  The  first  step  in  this  procedure  is  determination  of  the  val¬ 
ues  of  tangent  interface  stiffness  KJ  from  the  test  data  using  the  following 
equation: 
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(B2) 


K'st  = 


xf  —  x,— 1 


where 


X1'1 ,  xl  =  consecutive  shear  stress  readings 
A  A  1  =  consecutive  interface  displacement  readings 

£  5 


The  values  of  Kst  can  be  normalized  according  to  the  following 
expression: 


Ksn  = 


K's , 


\n> 


Yu-'  — 
\Pa 


(B3) 


where  Ksn  is  the  normalized  stiffness  of  the  interface.  The  parameter  tij  is 
the  stiffness  exponent  in  the  Clough  and  Duncan  (1971)  hyperbolic 
model,  which  is  determined  as  described  previously  in  this  appendix. 

Figure  B8  is  a  diagram  of  normalized  stiffness,  represented  in  logarith¬ 
mic  scale,  versus  stress  level.  The  data  points  correspond  to  the  initial  load¬ 
ing  tests  performed  on  the  interface  between  dense  Blacksburg  Sand  and 
concrete. 

G6mez,  Filz,  and  Ebeling  (2000)  propose  three  adjustments  to  the  hyper¬ 
bolic  parameter  values.  The  first  adjustment  consists  of  finding  the  value 
of  n-  that  minimizes  the  scatter  of  the  data  points  in  the  diagram  of  normal¬ 
ized  stiffness.  Ideally,  the  adjusted  value  of  would  produce  data  points  in 
the  diagram  that  lie  on  a  single  trendline  (Gomez,  Filz,  and  Ebeling  2000). 
In  practice,  the  optimal  value  of  is  that  which  minimizes  the  scatter  of 
the  data  points. 

For  the  interface  between  dense  Blacksburg  Sand  and  concrete,  it  was 
found  that  no  adjustment  was  required  for  the  nj  value  determined  pre¬ 
viously  in  this  appendix. 

The  second  and  third  adjustments  require  plotting  the  normalized  hyper¬ 
bolic  stiffness  Km,  which  is  calculated  from  the  following  equation: 

Ksn  =  Kr(\-RrSLj  (B4) 
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Figure  B8.  Adjustment  of  the  value  of  stiffness  number  Kt  for  the  dense 
Blacksburg  Sand-to-concrete  interface 

If  the  response  of  the  interface  is  hyperbolic,  accurate  values  of  Kj  and 
Rjj  would  produce  a  Ksn  versus  SL  diagram  that  fits  closely  the  normalized 
stiffness  diagram  of  the  test  data.  The  dashed  line  in  Figure  B8  corre¬ 
sponds  to  the  values  of  Ksn  before  adjustment  of  the  hyperbolic  parameter 
values  determined  previously  in  this  appendix.  It  is  seen  that,  although  it 
shows  good  agreement  with  the  test  data,  the  dashed  line  can  be  adjusted  to 
slightly  improve  the  accuracy  of  the  model.  By  reducing  the  value  of  Kj,  an 
excellent  agreement  with  the  test  data  is  achieved.  No  adjustment  was  nec¬ 
essary  for  the  value  of  R ^  determined  previously. 

Table  B1  summarizes  the  values  of  the  hyperbolic  parameters  for  the 
interface  between  dense  Blacksburg  Sand  and  concrete. 


Table  B1 

Hyperbolic  Parameter  Values  for  the  Interface  Between  Dense 

Blacksburg  Sand  and  Concrete 

Parameter1 

Parameter  Value 

K, 

23000  (adjusted) 

h _ 

0.8 

0.76 

8,  deg 

31.6 

1  Parameters  are  listed  and  defined  in  the  notation  (Appendix  C). 
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Appendix  C 
Notation 


B  Bulk  modulus  of  soil 

Cc  Coefficient  of  curvature 

Ck  Interface  stiffness  ratio 

Cu  Uniformity  coefficient 

Dr  Relative  density 

D !  Thickness  of  the  backfill  above  the  hydrostatic  water  table 

Thickness  of  the  submerged  backfill  above  the  heel  of  the 
wall 

Dio  30  60  Particle  size  diameter  corresponding  to  10,  30,  or  60  percent, 
respectively,  passing  in  the  grain  size  distribution  curve 

£.  Initial  tangent  (Young’s)  modulus  of  the  soil 

Gs  Specific  gravity 

I  Correction  factor  for  inclination  of  the  stress  path 

K  Modulus  number 

Kb  Bulk  modulus  number 

KAine  Line  joining  the  points  in  the  p'-q  plane  that  correspond  to 

failure 

Kj  Dimensionless  interface  stiffness  number  for  initial  loading 

Kq  At-rest  pressure  coefficient 

K  Initial  shear  stiffness  of  the  interface 

SI 
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ci 


Normalized  shear  stiffness  of  the  interface 


C2 


K, 


ts 


k  ys 

sn 


K 


St 


K.. 


K 

urj 


m 


ml 


"j 

P' 

Pa 

q 

Rf 

Rfj 

SL 

SL0 

SLts 

SVS 


a 


• max ’  ^ min 
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Transition  stiffness  number 
Yield  stiffness  number 

Interface  tangent  stiffness  for  vertical  stress  paths  (as 
defined  in  Clough  and  Duncan  1971  hyperbolic  model 
for  interfaces) 

Interface  tangent  stiffness  for  stress  paths  of  any  orientation 

Unload-reload  modulus 

Unload-reload  stiffness  number  for  interfaces 

Bulk  modulus  exponent 

Stiffness  degradation  parameter 

Modulus  exponent 

Interface  stiffness  exponent 

(cr/  +  o/)/2 

Atmospheric  pressure  =  101.3  kPa 
Shear  direction  parameter;  (Oj-o3)/ 2 
Failure  ratio  for  soils 
Failure  ratio  for  interfaces 
Stress  level 

Stress  level  at  the  origin  of  unloading-reloading 
Transition  stress  level 

Stress  level  for  current  position  of  yield  surface 
Scaling  factor  for  unloading-reloading 
Maximum  and  minimum  density,  respectively 
Unit  weight  of  water  (9.8  kN/m3) 

Peak  interface  friction  angle 
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A 


S 


A  U1  A  1 

“s  ’  ff 


A 


.so 


Ac() 


e 

ev 

e 

V 


v 

ynom 

°h 


0 


n 


°1 

°r 

°3 

°3> 

(07  -  o3) 

(Oj  -  Gj)f 

«*i  -  <* Ait 


Residual  interface  friction  angle 

Displacement  along  the  interface 

Consecutive  interface  displacement  readings 

Interface  displacement  at  the  origin  of  unloading-reloading 

Interface  displacement  to  peak 

Interface  displacement  to  residual 

Reduction  in  the  peak  secant  friction  angle  value  for  a 
10-fold  increase  in  03' 

Axial  strain 

Volumetric  strain 

Angle  between  the  stress  path  direction  and  the  x-axis 

Poisson’s  ratio 

Nominal  Poisson’s  ratio 

Horizontal  pressure 

Normal  stress  acting  on  the  interface 

Normal  stress  corresponding  to  point  TS  where  the  stress 
path  intersects  a  transition  surface 

Normal  stress  at  the  origin 

Major  principal  total  stress 

Major  principal  effective  stress 

Minor  principal  total  stress 

Minor  principal  effective  stress 

Deviator  stress 

Deviator  stress  at  failure 

Asymptotic  deviator  stress 


Appendix  C  Notation 


C3 


X 


Interface  shear  stress 


C4 


x1'1,  x' 


X 


ts 


xf 

Xo 

Xult 

<t> 

¥ 

^cv 

<t>* 


Consecutive  interface  shear  stress  readings 

Shear  stress  at  point  TS  where  the  stress  path  intersects  a 
transition  surface 

Interface  shear  strength 

Interface  shear  stress  at  the  origin  of  unloading-reloading 

Asymptotic  interface  shear  stress 

Peak  secant  internal  friction  angle  of  the  soil 

Effective  peak  secant  friction  angle  of  the  soil 

Friction  angle  at  a  strain  of  15  percent 

Peak  secant  friction  angle  at  a  confining  pressure  of 
101.3  kPa  (1  atm) 
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